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O : ABSTRACT 

CL|| We consider the simplest and most economic version among the proposed non-minimal 
QhI supersymmetric models, in which the /i-parameter is promoted to a singlet superfield, 
r-| , whose all self- couplings are absent from the renormalizable superpotential. Such a par- 
ticularly simple form of the renormalizable superpotential may be enforced by discrete 
^ ■ /^-symmetries which are extended to the gravity-induced non-renormalizable operators as 
' well. We show explicitly that within the supergravity-mediated supersymmetry-breaking 
scenario, the potentially dangerous divergent tadpoles associated with the presence of the 
gauge singlet first appear at loop levels higher than 5 and therefore do not destabilize the 
gauge hierarchy. The model provides a natural explanation for the origin of the /i-term, 
without suffering from the visible axion or the cosmological domain-wall problem. Focus- 
ing on the Higgs sector of this minimal non-minimal supersymmetric standard model, we 
calculate its effective Higgs potential by integrating out the dominant quantum effects due 
to stop squarks. We then discuss the phenomenological implications of the Higgs scalars 
predicted by the theory for the present and future high-energy colliders. In particular, we 
find that our new minimal non-minimal supersymmetric model can naturally accommodate 
a relatively light charged Higgs boson, with a mass close to the present experimental lower 
bound. 
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1 Introduction 



In the well-established Standard Model (SM), the generation of gauge-invariant, renormal- 
izable masses for the observable fermions, e.g. the electron and the t quark, and for the W 
and Z bosons is achieved through the so-called Higgs mechanism. Most interestingly, the 
Higgs mechanism itself predicts inevitably the existence of a fundamental scalar, known as 
the Higgs boson. Recently, experiments at LEP2 have intensified their searches for directly 
observing the yet-elusive Higgs boson. Their latest analyses show that its mass must be 
larger than 113.3 GeV at the 95% confidence level (CL) |[l[]. At the same time, electroweak 
precision data place an upper bound of the order of 240 GeV on the Higgs-boson mass f^. 

So far, we have no much evidence to suggest that the underlying structure of the Higgs 
potential is indeed that of the SM or it already contains components of a more fundamen- 
tal theory which is about to be unraveled in the next-round experiments. In particular, it 
is known that the SM cannot adequately address the problem of gauge hierarchy, which 
is related to the perturbative stability of radiative effects between the electroweak scale 
and the Planck or grand unification scale. An appealing solution to this problem may be 
achieved by means of supersymmetry (SUSY). In order that SUSY theories avoid reintro- 
ducing the problem of gauge hierarchy, they must be softly broken at a relatively low scale 
MsvsY ~ rrit of the order of 1 TeV, in agreement with experimental observations. 

The minimal supersymmetric extension of the SM, also called the Minimal Supersym- 
metric Standard Model (MSSM), predicts a very constrained two-Higgs-doublet potential 
at the tree level, whose quartic couplings are determined by the well-measured SU(2)i and 
U(l)y gauge couplings gw and g'. As a consequence, the lightest neutral Higgs boson is 
always lighter than the Z boson at the tree level. Nevertheless, radiative corrections to 
the effective Higgs potential are significant and extend the above mass upper bound to 
110 (130) GeV for small (large) values of the ratio of Higgs vacuum expectation values 
(VEV's) tan/3 ^ 2 (> 15) 0. Thus, a large portion of the parameter space of the MSSM 
has been already excluded by the current LEP2 experiments at CERN. Moreover, the up- 
graded Tevatron collider at Fermilab will have a much higher reach in discovering heavier 
Higgs bosons with SM-type couplings and masses up to 140 GeV and therefore will provide 
a unique test for the viability of the MSSM. 

On the basis of the above strong experimental bounds on the lightest Higgs-boson 
mass in the MSSM (especially for low values of tan/?), it would be rather premature to 
infer that realizations of low-energy SUSY in nature have a rather limited range. In order 
to reach a more definite conclusion, it is very important to further analyze the Higgs sectors 
of minimally extended scenarios of the MSSM. An additional reason for going beyond the 
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MSSM is the so called /i-problem. The superpotential of the MSSM contains a bilinear 
term —HH1H2 involving the two Higgs-doublet superfields Hi and H2, known as the fi- 
term. Although /i is naturally of the order of the Planck scale Mp, it is actually required 
to be many orders of magnitude smaller of order Msusy for a successful Higgs mechanism 
at the electroweak scale. Many scenarios have been proposed in the existing literature to 
account for the origin of the /x-term, albeit all in extended settings [Q. 

A simple SUSY extension of the MSSM, which one might have thought of considering 
to address the /i-problem, would be to elevate the /i-parameter to a dynamical variable by 
means of a gauge-singlet chiral superfield S, couple the latter to Hi and H2 as XSH1H2 
and arrange that S somehow develops a VEV of the order of Msusy ~ ""^t- However, 
this minimally extended scenario possesses a global U(l) Peccei-Quinn (PQ) symmetry, 
whose spontaneous breakdown gives rise to a phenomenologically excluded axion. The 
most popular way in the literature of removing the unwanted PQ symmetry is to break 
the latter explicitly by adding the cubic self-coupling to the superpotential. The 

resulting model has been termed the Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) p. Unfortunately, the NMSSM is also plagued by its own problems. The cubic 
self-coupling of S leaves invariant a subgroup of U(1)pq, namely the discrete Z3 symmetry, 
whose subsequent spontaneous breakdown gives rise to the formation of cosmologically 
catastrophic weak-scale domain walls [^,0. 

Another well-known problem that a model of low-energy physics involving light gauge- 
singlets has to face is the destabilization of the gauge hierarchy through the generation of 
at least quadratically divergent tadpoles for the singlet [Q. In the context of = 1 
supergravity, which is spontaneously broken by a set of hidden sector superfields, even if one 
assumes no other scale between Msusy and Mp, the simple presence of gravity- induced non- 
renormalizable operators in the superpotential and the Kahler potential is able to generate 
such tadpoles Using the Planck mass Mp as a physical cut-off energy, such divergences 
contribute tadpole terms of order (l/167r^)"MpM|usY'S' to the effective potential, where 
n indicates the loop level at which the tadpole divergence appears. It is obvious that for 
small values of n, e.g. n < 4, the generated tadpole terms lead generically to unacceptably 
large values for the VEV of S (the scalar component of S), thereby destabilizing the gauge 
hierarchy. 

In the case of the aforementioned extensions of the MSSM, the problem of destabiliza- 
tion does not occur as long as the U(1)pq or Z3 symmetries are imposed on the complete 
set of non-renormalizable operators as well. However, any attempt to break these unwanted 
symmetries through a subset of non-renormalizable operators would, as an immediate con- 
sequence, destabilize the weak scale. This aspect has been emphasized in Ref. [0,0], in 
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connection with the ^3 symmetry of the NMSSM. 

Recently, it has been reahzed that the unwanted ^3 and U(1)pq symmetries present in 
the corresponding supersymmetric extensions of the MSSM could be effectively broken not 
by the non-renormalizable operators themselves, but rather by the tadpoles generated by 
them ||ll],|l2|. For the 23-symmetric extension of the MSSM, harmless tadpole terms of order 
(l/167r^)"'M|ugY5', with 2 < n < A, were sufficient for the breaking of the Z3 symmetry. 
For the PQ-symmetric extension instead, it was necessary that the harmful tadpoles of 
order (l/167r^)"MpM|ugY5' (using Mp as a cut-off scale) be generated at a sufficiently high 
loop level n, with 5 < < 8. In the Z3 case, the harmful tadpoles were forbidden by 
imposing on the operators of the non-renormalizable superpotential and Kahler potential 
the Z2 i?-symmetry of the cubic superpotential, under which all superfields as well as the 
superpotential flip sign. However, the desirable form of the renormalizable superpotential 
was enforced by imposing a larger group, namely the product of the Z2 matter parity with 
a Z4 i?-symmetry [0. In the U(1)pq case, a Z^ i?-symmetry proved sufficient to enforce 
the desirable renormalizable superpotential and postpone the appearance of the harmful 
divergent tadpoles until the sixth loop order [|12|. Thus, in both cases the breaking of the 
unwanted symmetries was successfully implemented without jeopartizing the stability of 
the electroweak scale and without generating new cosmological problems. 

In the present paper, we shall study in detail the new minimal supersymmetric ex- 
tension of the MSSM, in which the linear, quadratic and cubic terms involving the singlet 
superfield S itself are absent from the renormalizable part of the superpotential. Hereafter, 
we shall call such a supersymmetric extension the Minimal Non-minimal Supersymmetric 
Standard Model (MNSSM). In particular, we shall explicitly show that with the imposition 
of the discrete Z^ and Zi /^-symmetries on the complete superpotential and on the Kahler 
potential of the corresponding supergravity models, the potentially dangerous tadpole di- 
vergences first appear at the six- and seven- loop levels, respectively, and hence are naturally 
suppressed to the order of M|usy S. Evidently, the resulting model constitutes the simplest 
and most economic version among the non-minimal supersymmetric models proposed in 
the literature. In order to properly study the properties of the Higgs bosons predicted by 
the theory, we will calculate the effective Higgs potential by taking into account the domi- 
nant stop-loop effects. Finally, we shall analyze the phenomenological implications of the 
MNSSM for direct Higgs-boson searches at the LEP2 and the upgraded Tevatron colliders. 

The organization of the paper is as follows: in Section 2 we describe the Higgs sector 
of the MNSSM and show that harmful tadpole divergences first appear at the six- and 
seven- loop levels, after the aforementioned discrete and Z-j /^-symmetries are respec- 
tively imposed on the theory. Technical details of the argument are relegated to Appendix 
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A. In Section 3 we compute the effective Higgs potential by integrating out the domi- 
nant radiative effects due to stop squarks, from which we derive the CP-even and CP-odd 
Higgs-boson mass matrices. In Section 4 we investigate the theoretical differences of the 
Higgs-boson mass spectrum between the MNSSM under consideration and the frequently- 
discussed NMSSM. In Section 5 we present numerical estimates of the Higgs-boson masses 
and their couplings associated with the Z boson in these two models, and discuss the phe- 
nomenological implications of the MNSSM Higgs sector for the direct Higgs-boson searches 
at LEP2 and for the upcoming searches at the upgraded Tevatron collider. Section 6 
contains our conclusions. 



2 MNSSM: Symmetries and stability of 
the electroweak scale 

In this section we shall consider the simplest extension of the MSSM, the MNSSM, within 
the context of = 1 supergravity spontaneously broken by a set of hidden sector fields at 
an intermediate scale. In the MNSSM, the /i-parameter is promoted to a dynamical chiral 
superfield,Q with the linear, quadratic and cubic terms involving only the singlet superfield 
S being absent from the renormalizable superpotential. Such a particularly simple form of 
the superpotential may be enforced by discrete i?-symmetries, e.g. and Z^^ which are 
extended to the non-renormalizable parts of the superpotential and the Kahler potential as 
well. Adopting the standard power counting rules QJl^l, we shall show that in such A^ = 1 
supergravity scenarios, the potentially dangerous tadpole divergences are suppressed by 
loop factors l/(167r^)" of order = 6 and higher, and therefore do not destabilize the 
gauge hierarchy. Technical details are given in Appendix A. 

The renormalizable superpotential of the MNSSM under discussion is given by 

W^ren = klH^lT^LE + hdH^lT2QD + KQ^IT2H2U + XSHJ1T2H2. (2.1) 

where T2 is the usual 2x2 Pauli matrix. In Eq. (|2.1|), the Higgs superfields. Hi and H2, 



as well as the quark and lepton chiral multiplets, Q and L, are SU(2)i-doublets, while the 
remaining superfields U, D and E are singlets under SU(2)i. The chiral multiplets also 
carry the following hypercharges: 

U(l)^ : Hi (-1), (1), ^ (0), Q (1/3), U (-4/3), D (2/3), L (-1), E (2) , (2.2) 



*An earlier suggestion along these lines was discussed in |T3[ 
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where the hypercharge of each superfield is indicated within the parentheses. In addition 
to the baryon (B) and lepton (L) numbers, the renormahzable superpotential PVren respects 
the global U(l) PQ and R symmetries: 

U(l)pQ : ifi(l), i?2(l), ^(-2), g(-l), [/(O), D(0), i(0); 
U(l)« : ifi(O), i?2(0), 5(2), Q(l), f/(l), W^rcn (2) • (2.3) 

Note that IVrcn has charge 2 under U{l)j^. The symmetry group U{l)j^ is non-anomalous 
with respect to QCD interactions, but gets broken by the soft SUSY-breaking trilinear 
couplings down to its maximal non-i? Z2 subgroup which becomes the known matter- 
parity. Instead, the anomalous symmetry U(l)pQ remains unbroken by the soft SUSY- 
breaking terms. Neglecting QCD-instanton effects, U(l)pQ will remain unbroken, unless a 
gravity- induced tadpole operator linear in 5* gets generated from the non-renormalizable 
sector of the theory. The tadpole operator generically contributes to the effective potential 
a term 

^tad ~ p^MpMsV^ + h.c., (2.4) 

where n is the loop level at which the tadpole divergence occurs, using the Planck mass Mp 
as an energy cut-off. The tadpole term Vtad together with the soft SUSY-breaking mass 
term MI^syS*S lead to a VEV for the singlet field S of order (l/167r2)"Mp. To avoid 
destabilizing the gauge hierarchy, one must require that {S) ~ Msusy ~ (l/167r^)"'Mp, 



with Msusy ~ 1 TeV. This requirement can only be fulfilled for sufficiently high values of n, 
i.e. for n > 5. Finally, we should remark that the full renormahzable Lagrangian, including 
the tadpole term, preserves the B and L numbers. However, the quantum numbers B and L 
may be violated by certain non-renormalizable operators, which are hopefully of sufficiently 
high order in order not to upset the laboratory limits on proton instability. We can therefore 
conclude that the renormahzable superpotential H^ren of Eq. ( p.l|) supplemented with a 
sufficiently suppressed tadpole for the singlet S leads to a model without any obvious 
phenomenological or cosmological problem. 

One may now wonder whether there exists a symmetry giving rise to the above- 
described model that includes a tadpole term for S of the desirable order. To address this 
question, let us consider the global symmetry defined as a linear combination R' = 3/2+ PQ 
of U(1)r and U(1)pq, with 

U(l)«, : ifi(l), H2{1), S(4), Q(2), f/(3), D{3), L{2), E{3), W,en (6) . (2.5) 

Observe that the imposition of U(l)^, is sufficient to ensure the form ( p.l| ) for Wren- We 
should now examine whether U(l)^, also allows the generation of a tadpole term. The 
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(a) (b) 

Figure 1: Typical harmful tadpole divergences at the (a) six- and (b) seven- loop levels. 

symmetry group U(l)^, is explicitly broken by the trilinear soft SUSY-breaking interactions 
down to its maximal non-R subgroup Zq which is isomorph (equivalent) to the product 
group Z2 X Z^. The symmetry Z2 is essentially the ordinary matter-parity, under which 
the tadpole remains invariant. Instead, the symmetry Z^ is broken by the tadpole of S. 
Consequently, a tadpole term can only be generated if the whole symmetry group U(l)ij' 
or one of its subgroups that includes Z3 is violated by the higher-order non-renormalizable 
operators. 

The above arguments seem to suggest that the symmetry we are looking for is likely 
to be a subgroup of U(l)/j/ which is sufficiently large to enforce the form of W^-en given by 
Eq. (|2.1| ), but does not contain the Z3 subgroup of U(l)j?/. Subgroups of \J{1)r' obeying 
the above criteria are the discrete /^-symmetries Z^ and Z^^. 

Let us first consider the Z^ case. Under Z^, the chiral multiplets as well as the 
superpotential PVren transform as follows: 

Z^ : (H,, %) ^ u: {H,, %) , (Q, L) ^ (Q, L) , (f/, D, E) ^ uj^ {U, D, E) , 

S^UJ^S, IVren^^W^ren, (2.6) 

with uo = exp(27ri/5) and = 1. The discrete R -symmetry Z^ is imposed on the complete 
superpotential and Kahler potential. By means of standard power counting rules governing 
the harmful tadpole divergences, it can be shown that harmful tadpoles first appear at the 
six-loop level. As can be seen from Fig. |l](a), a typical harmful six-loop tadpole diagram 
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can be induced by appropriately combining the non-renormalizable operators of the Kahler 
potential 



K2 = 1^2 + h.c, K5 = K5 + h.c, (2.7) 

and four times the renormalizable term X SH^iT2H2 of the superpotential ( p.l|) . The an- 
alytic steps of the argument are presented in Appendix A. Thus, the induced harmful 
divergent tadpole term has the form 

^tad ~ l^g^MpMlusY^ + h.c. (2.8) 

From Eq. (|2.^ ), it is easy to see that the tadpole term is of order (1 TeV)^, e.g. for 
/t2 ~ K5 ~ 0.1, A ~ 0.6 and Mgusy ~ 1 TeV, and does not destabilize the gauge hierarchy. 

In the case, the harmful tadpole divergence occurs at one loop-order higher, 
namely at the seven-loop level, so the generated tadpole terms can naturally be as low 
as (100 GeV)^. In detail, under this new discrete R -symmetry, the superfields and Wren 
transform in the following way: 

: {Hi, %) -^00 {Hi, %) , (Q, L) ^ uj^ (Q, L) , (f/, D, E) ^ uo"" {U , D, E) , 

S^UJ^S, Wren^w'Wren, (2.9) 

with u = exp(27rz/7) and cj'' = 1. Following the same line of steps as above, we impose the 
discrete /^-symmetry on the complete superpotential and Kahler potential. Based on 
standard power counting rules, we show in Appendix A that the potentially harmful tadpole 
divergences first appear at the seven-loop level. A typical harmful tadpole diagram at seven 
loops is displayed in Fig. |l](b), and can be obtained by combining the non-renormalizable 
operators of the Kahler potential 

and four times the renormalizable term A SHfiT2H2 of Wren- The size of the so- generated 
tadpole term may be estimated as 

^tad ~ -^^^MpMi^^yS + h.c ~ (ITeV)xMsV^ + h.c, (2.11) 

for ~ Kg ~ 1 S'lid A ~ 0.6. If K^l^ ~ Kg ~ 0.1 and Mgusy ~ 1 TeV, the size of V^ad can 



be as low as (0.2 TeV) 
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We conclude this section by noticing that although the discrete i?-symmetries 
and Zl^ do not contain the usual Z2 matter parity, they still prohibit the presence of all 
dimension ci = 4 B- and L- violating operators as well as the dangerous B- and L- violating 
operators QQQL and UUDE of dimension 5. However, the symmetries Z^ and 2^ allow 
the L-violating operator LLH2H2 of c? = 5, which is able to generate Majorana masses 
for the light left-handed neutrinos. Moreover, Z^ allows the d = 5 L-violating operators 
SSLH2, SLLE and SLQD, whereas Z^ allows the d = 5 B-violating operator SUDD. 
Although these last operators are unable to lead by themselves to an observable proton 
decay, they still render the lightest supersymmetric particle (LSP) unstable. However, 
estimates based on naive dimensional analysis show that the LSP is very long lived with a 
lifetime larger than the age of the Universe and therefore safely qualifies to be a dark-matter 
candidate. Of course, the LSP can be made absolutely stable by the additional imposition 
of the Z2 matter parity. 

3 The Higgs sector of the MNSSM 

In this section we shall study the low-energy Higgs sector of the MNSSM. After discussing 
its tree-level structure, we will then calculate the one-loop effective Higgs potential by 
integrating out the dominant loop effects due to stop/top quarks, from which we derive 
analytic expressions for the Higgs-boson masses and their respective mixing angles. We 
shall then focus on the gaugino-Higgsino sector of the MNSSM, and briefly discuss possible 
laboratory limits on the would-be //-parameter due to the presence of a light quasi-singlet 
neutralino state. Finally, for our forthcoming phenomenological discussion in Section 5, we 
shall present the effective Higgs-boson couplings to the W and Z bosons. 

3.1 Higgs-boson masses at the tree level 

In addition to terms proportional to S, another effect of the tadpole supergraphs of Fig. 
|T] is the generation of terms proportional to Fs, namely to the auxiliary scalar component 
of S. As a consequence, the effective renormalizable Higgs superpotential of the MNSSM 
reads 

W^ffiggs = ^SH^^r2H2 + ^fMIusy^, (3.1) 

where C,f is a model-dependent constant. Moreover, the Lagrangian describing the soft 
SUSY-breaking Higgs sector is given by 

- Aoft = (e^Mlusv^ + h.c.) + m2$l$i + ml^l<!>2 + mlS*S 



9 



+ [XAx S^\tT2^2 + h.c.) , (3.2) 

where $i = 2x2$* and $2 are the physical bosonic degrees of freedom of Hi and H2, 
respectively. After including the relevant F- and D-term contributions in addition to the 
soft SUSY-breaking terms, we obtain the complete renormalizable Higgs potential of the 
model of interest 



C^y = {tsS + h.c.) + ml^l^i + m^$^$2 + mlS*S + (m^2$I<l>2 + h.c. 

+ (AAa^$I$2 + h.c.) - Ai(<l'I$i)2 - A2($5$2)' - A3($I$l)($5$2) 

- (A4 - A^) ($l$2)(<l>5$i) + A^ 5*5(<l>I$i + $^$2) , (3.3) 



with 



Ai = A2 = A3 = A4 = (3.4) 

Here, {g') is the coupling constant of the gauge group S'U{2)l (U(l)y). As was discussed 
in the previous section, the tadpole prefactor ^5 in Eq. ( |3.4|) is of order unity. However, 
the size of $,f crucially depends on the VEVs of the scalar components of the hidden-sector 
superfields that break SUSY Q . The VEVs remain unconstrained by the requirement that 
the breaking of SUSY takes place at some intermediate scale in the hidden sector, in which 
the F-terms of the respective hidden-sector superfields are involved. In case that some of 
the hidden-sector fields acquire VEVs of order Mp, the tadpole prefactors \^p\ and l^^l 
could be comparable. Otherwise, it is l^^l <^ \C,s\. In the following, we shall treat the ratio 
/ 1^5 1 as a free parameter which is always less than unity. 

We shall now derive the minimization conditions of the Higgs potential in Eq. ( p.3|) . 
Throughout the paper, we shall assume that CP is a good symmetry of the theory. Under 
this assumption, we can perform the following linear expansions of the Higgs fields about 
their VEVs: 




(3.5) 



The minimization conditions are then determined by the vanishing of the tadpole param- 
eters 



ml + (^-^ XAxVs + ml^ ) tp - Xivf 
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I (A3 + A4 - A^) vl + I \\l 



2 

dCv 



V2 



TTlr, 



XAxVs + m 



2 

dCv 



i\3 + K-\^)vl + -A^t;^ 



-Vs [mg + XA) 



V1V2 



V2 



Vs 



A2f2 



Vs 



(3.6) 

(3.7) 
(3.8) 



with V 



vf + v^ 



— 2M]v/g^ and tf^ = V2/V1. Our earlier assumption of CP invariance 
entails that all kinematic parameters involved, e.g. A and Ax, are real, namely there are no 
explicit sources of CP violation in the theory. Also, it is important to remark that based 
on Romao's no-go theorem |T^, CP invariance cannot be broken spontaneously at the tree 
level in the MNSSM.Q 

It proves now convenient to perform a change of the weak basis for the charged and 
CP-odd scalars: 

— ^0 






(3.9) 



where S/? = V2/V and cp = vi/v, such that H'^ becomes the mass eigenstate of the charged 
Higgs boson, and and G° are the would-be Goldstone bosons which constitute the 
longitudinal degrees of freedom of the and Z bosons, respectively. 

Let us first consider the charged Higgs sector. In the newly defined weak basis of Eq. 
(^.9[), the tree-level mass of the charged Higgs boson may easily be computed by 

where ^ 

= --^Xvs (3.11) 

is the would-be /i-parameter of the MSSM. Here and in the following, we adhere the super- 
script (0) to a specific kinematic quantity in order to emphasize its tree-level origin, e.g. 



m 



12 



+ M, 



w 



^ \2 2 

- X V . 



(3.10) 



2(0) 



H+ 



Since the would-be Goldstone boson does not mix with other fields, the tree-level 
CP-odd mass matrix takes on the simple form in the reduced weak basis {a, as}'- 



M 



2(0) 



V 
Vs 



2(0) 



+ m 



12 



2(0) 



+ m 



12 



V 



+ 



\ts 



(3.12) 



^We find that this property persists, even if CP-conserving radiative effects mediated by large stop 
mixing are included in our model. 
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with 

^2(0) ^ ^2(0) _ + 1 . (3.13) 

In deriving the above form of Mp '''' , we have also considered the tadpole constraints given 
by Eqs. (pl6| ) -( ^^) . In the MSSM limit, which is obtained for vs ~ —\/2ts /m?s ^ with 
the would-be /i-parameter being kept fixed (A — 0), the mass eigenvalues of the CP-odd 
mass matrix can easily be approximated by 

mT ^ ^ — > (3-14) 

with Xts/fx > M^^^\ Furthermore, in the limit, in which the tadpole parameters 

ts vanish, the CP-odd mass matrix contains a massless state, i.e. a PQ axion, as a result 

of the spontaneous breakdown of the symmetry group f/(l)pQ. 

Taking into account the tadpole constraints of Eqs. (|3l6| )-(|3l8|), the tree- level CP-even 
mass matrix may be expressed in the weak basis {0i, 02, 0s} as follows: 

(Mf))n = cjMl + 4m5°), 

iMT% = {Mf% = -s,c,{Ml^'^ + Ml - AV), 



,2 



(Mf ))33 = ^.,c,(.,c,m5°) + ml,) + ^, (3.15) 

Vs ^ 



with Mz = ^gl + g'^v/2. In the MSSM limit, in which vs ~ -V2ts/ml > v with 
fixed, the Higgs-singlet components decouple from the tree-level CP-even mass matrix 
mI^^\ In this case, the heaviest Higgs boson Hs is predominantly singlet and has a squared 
mass M]j^^ ^ Xtsl jji] becomes mass degenerate with A, (cf. Eq. ( ^.14| )). 

Apart from the MSSM limit mentioned above, there exists a novel non-trivial decou- 
pling limit for the heavy Higgs sector in the MNSSM. This decoupling limit is obtained for 
large values of the tadpole parameter \ts\-, where all other kinematic parameters are kept 
fixed. In this case, the Higgs states A, and are singlets, i.e. A, = as and = 0^, 
and so decouple from the remaining Higgs sector while being degenerate in mass, i.e. 
M^^^"'* ^ M^^"* ^ Xts/fJ'- An immediate consequence of this is the relation 

Mlf ^ M^W, (3.16) 
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where M^^^^ is defined in Eq. ( p.l3| ). Most importantly, in this hmit the structure of the 
low-energy Higgs sector, although reminiscent of, is not identical to that of the MSSM. For 
example, as opposed to the MSSM limit, the terms proportional to A^f ^, which occur in the 
CP-odd and CP-even mass matrices of Eqs. ( |3.12| ) and ( p.l5| ), do not necessarily vanish in 
the decoupling limit due to a large tadpole. Thus, contrary to the MSSM, Eq. (|3.16| ) implies 
that for large values of A, e.g. A ~ g^^, the charged Higgs-boson mass can become 

even smaller than the mass Mi°^ of the non-decoupled CP-odd scalar. As we will see in 
Section 5, this last fact plays a very important role in lowering the mass of the MNSSM 
charged Higgs boson up to its experimental lower bound, i.e. up to Mh+ ~ 80 GeV |[T|Jl5 



Moreover, in Section 4 we shall see that this new non-trivial decoupling limit due to a large 
tadpole parameter Ats//i is only attainable in the MNSSM, and no analogue of this exists 
in the NMSSM. 



As in the MSSM |T6|, an upper bound on the mass of the lightest CP-even Higgs 
boson in the MNSSM with large \ts\ may easily be derived in the decoupling limit of a 
heavy charged Higgs boson, i.e. for 

^ » Mj(? » M|. (3.17) 

In this limit, in addition to A2 and H^, the Higgs scalars Ai and H2 decouple from the 
lightest Higgs sector as well, and are almost mass degenerate with the charged Higgs boson 
H^. After taking into consideration the heavy //^-decoupling limit of Eq. (|3.17|) , the mass 
of the lightest CP-even Higgs state Hi is found to satisfy the inequality 

M^°) < M| ( cos2 2/5 + sm'2p) . (3.18) 

Note that as opposed to the MSSM where A = 0, can be larger than Mz in the 

MNSSM, especially for small values of tan /3. This prediction is very similar to the one 
obtained in the frequently-discussed NMSSM [^. However, as tan/5 increases, e.g. for 
tan /5 ^ 5, the A-dependent term in Eq. ( |3.18| ) becomes negligible. Thus, in the large-tan /3 



case, the upper bound on the f/i-boson mass is almost identical to the one obtained in 
the MSSM. Finally, an important property of the MNSSM is that the tree-level neutral 
Higgs-boson masses satisfy the equality 

Mf^ + Mf^ + Mf^ = Ml + mJ°) + Mlf . (3.19) 

It is interesting to notice the striking similarity of the above mass sum rule with the 
corresponding one in the MSSM [|l6l, in which case the mass terms mJ°^ and M^f are not 



present in Eq. ( [3.191 ). The above observation allows us to advocate that the structure of 
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the MNSSM Higgs sector departs indeed minimally from that of the MSSM. Nevertheless, 
exactly as happens in the MSSM 0], the tree- level Higgs sector receives sizeable quantum 
corrections due to stop squarks, leading to a violation of the mass sum rule 

3.2 One-loop effective potential 

We shall now calculate the dominant one-loop corrections to the effective potential due to 
top (t) and scalar-top (t) quarks. As a good approximation, we neglect the one-loop D- 
term contributions as well as bottom (6) and scalar-bottom (6) quark effects by assuming a 
vanishing 6-quark Yukawa coupling, i.e. = 0. The above approximations are reasonable 
for relatively small values of tan/3, e.g. tan/? < 10, where the MNSSM predictions for the 
lightest Higgs sector are expected to deviate considerably from the ones obtained in the 
MSSM. 

The interaction Lagrangians relevant for the computation of the one-loop effective 
potential are given by 

-Acrmion = htQLir2^tji + h.C, 

-Cf = h^tl^ltr^QLl^ + {XhtSQitT^^tn + h.c.) + , 
-Csoft = M^QIQl + M^i^^tR + [htAtQyT2<^liR + h.c), (3.20) 

where Ql = {ti, bi)'^ and Ql = {ti, &l)^ are the bosonic and fermionic degrees of freedom 
of the third-generation left-handed quark superfield. 

Equipped with the Lagrangians in Eq. (|3.20| ), we can now derive the Higgs-dependent 
t and i masses. Thus, the squared t-quark mass in the Higgs background is given by 

ml = hi ^1^2 ■ (3.21) 

The corresponding background-dependent stop masses may be determined from the 3x3 
squark mass matrix, which is expressed in the weak basis {Ql, tn] as: 

~2 _ I M^l2 + /l? ($5$2l2 - $2$2) htAtlT2<i>*2 + Xht SiT2<i>l 

The squared squark-mass matrix Ad"^ has three mass eigenvalues. For = 0) these are 
given by the two squared Higgs-dependent t-quark masses 

= ^ (^Q + + 2/^? + (-) ^/iM^Q - + ^h'MtC^^ + XS*^l\^ ) 

(3.23) 
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and by the squared left-handed sbottom mass m?^ 



Mq, where 0°2 = ("^1,2 + 0i,2 + 



mi^2)/"\/2 are the neutral parts of $1,2- 

In the MS scheme, the one-loop Coleman-Weinberg effective potential [jl^ may be 
expressed in terms of the relevant squared Higgs- dependent masses m|, m|^, rh^^ and mf^ 
as follows: 



V 



327r2 



(3.24) 



k=ti,t2,bL 

where —Cy is the bare Higgs potential given by Eq. With the help of Cy, we can 

now compute the radiatively corrected mass of the charged Higgs boson by means of the 
relation 

1 / d^Cy 



2(0) 



H+ 



+ AM 



H+ 



M 



2(0) 
H+ 



E 

k=ti,t2,bL 



TJlr 



2 



(3.25) 



where M'^+' is the tree-level contribution and (ml) = m| . Following the procedure outlined 
in [|18| ,0 we find 



+ 



2 
ti 



+ 



m 



t2 



2 (m? . - 

^ tl(t2) 



it2) "^L)H-<^ 



2 (m? — m? ) (m? — ml ] 



(3.26) 



Then, the one-loop correction to Mjj+ , AM^+ , is given by 



In 



2 2 



+ 



(3.27) 



with 



g{m\,ml) 



2 



2. 



(3.28) 



In Eq. (|3.27] ), the quantity 5rem summarizes the remaining Q^-independent corrections: 



3^t;2 
647r2 



2 



In 



t2 



2 



t2 



+ 



.) V 

2 
t2 



m 



ti 

^ht _ 

327r2 m? + m? 

ii 12 



bL 



2 2 
— m= 

1-2 bL 



bL 

In 



2 



i2 



■l-A similar procedure was also followed in |19[ 



(3.29) 
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As we will see below, the presence of Sj-em gives rise to a modification of the tree-level relation 
between M|^+ and in Eq. ( 3.13|) . Nevertheless, it can be estimated from Eq. (|3.29 ) 



that this modification, which scales quadratically with the /i-parameter, is insignificant for 
almost all relevant values of fi of interest to us, i.e. for l/x/m^-J < 2. 

We now calculate the one-loop radiative shift AMp to the CP-odd Higgs-boson mass 
matrix. The analytic result may be completely expressed in terms of AM^ = AM^+ — 5rem 
as 

1 ^ SpC^ 

AMI = AM^ I ^ v' 2 2 \ ■ (3-30) 



It is easy to see that the one-loop radiative shift may be entirely absorbed into the tree- 
level mass matrix in Eq. (|3.12| ), after performing an one-loop re-definition of M^^^\ namely 
= M^^^-* + AM^. After this re-definition, the tree- level mass relation in Eq. (|3.13| ) gets 
radiatively corrected as follows: 

= Ml+ - Mir + ^AV - (5rem. (3.31) 

The one-loop Born-improved CP-odd mass matrix, M|, = Mp^°^ + AM|,, may be 
diagonalized through an orthogonal transformation of the weak fields 

, with = i '^"^^ \ . (3.32) 

\ — sin6'^ cos6'a / 

The CP-odd fields Ai and A2 are the mass eigenstates of M|,, with squared masses 





= i (TrM|, - (+) ^Tt^MI - 4detM|) . (3.33) 



The mixing angle 9a relating the weak to the mass eigenstates is uniquely determined by 

\{Mi),,\ . mi),, - Mij 

cost^4 = = , smt^/ — 



(3.34) 

Finally, we calculate the radiative corrections AM| to the CP-even Higgs-boson mass 
matrix. The individual matrix elements of AM| are given by 



t2' 

{AMD,, = (AM|)2i = -wAMf - 



167r2 



2 2 \ 2 
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Vs 
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(3.35) 



where Xt = At — n/tp. Again, we find that almost the entire Q^.j^ependence of the 
radiatively-corrected CP-even mass matrix given in Eq. ( p.35|) can be absorbed into M^ 
by an one-loop re-definition of M^^^\ An exception to this is the mass-matrix element 
(AM|)i3. The Q^-dependence of the {13} element can be eliminated by the $2-"wave- 
function counter term (CT) which is contained in the A-parameter. 

To make this last point explicit, we shall apply the non-renormalization theorem of 
the superpotential to the coupling \SHjiT2H2 in Eq. (p.l|). Since this operator does not 
receive any ultra-violet (UV) infinite radiative corrections to all orders, the wave-functions 
of S, Hi and H2, denoted as Z-^, Zcj and Z-Tj must cancel against the CT of A, 5X, that is 



5\ 



iz-J 

V 5 



/22— i/a^— 1/2 



H2 



l)A = -\{6Z^ + 6Z^^+ 8Z^^)\ 



(3.36) 



where Z}J'' = 1 + \5Zz, with z = S, Hi, and H2- Since only the wave-function of H2 
receives quantum corrections due to top quarks, Eq. ( |3.36D becomes 

SXhl 



--6Z- A 
2 ^2 



327r2 



^2 fmf 
In — ^ 



(3.37) 



Here, we have implicitly assumed that the coupling A is renormalized at the scale Q' 



ml Returning now to the bare Higgs potential in Eq. ( p.3| ), we see that the operator 
A^ (S*^)($l$i) induces the CT 2XSX, which gives rise to a corresponding CT in the tree- 
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level mass-matrix element (Mi 



{SM; 



I)l3 



V vs 



fx 



3hf / vcp 
8^ \ 



mt 



(3.38) 



Adding the CT ((5M|)i3 to the one-loop result (AM|)i3, we readily see that gets sub- 
stituted by m^. Finally, it is not difficult to convince ourselves that there are no analogous 
5A-dependent CTs for the operators {S*S){^1^2) and ($|$i)($5$2), as they are ex- 
actly canceled by the wave-function renormalization of $2- 

The one- loop radiatively corrected CP-even mass matrix, M| = M'^!"^^ + ^M"^, is 
i of a 3 X 3 orthogonal matrix , i.e. 

(0^)^M|0^ = diag(M^^, M^^, 

Mfj,^. Under this 
related to the mass eigenstates through 



diagonalized by means of a 3 x 3 orthogonal matrix O , i.e. 



(3.39) 



with Mfj_^ < Mfj^ < Mfj,^. Under this orthogonal transformation, the weak states are 



Ui ] 




( 


02 


= 0^ 




V 05 ) 







(3.40) 



The entries of can be calculated analytically by solving the third-order characteristic 
equation of M|. The procedure of deriving analytic expressions for the elements of is 
very similar to the one presented in Appendix B of [^^, and we will not repeat it here. 



3.3 The Higgsino sector 



In addition to the Higgs sector, the Higgsino (or neutralino) sector of the MSSM gets 
minimally extended in the MNSSM due to the presence of the neutral SUSY partner of the 
complex scalar singlet the singlino s. Instead, the tree-level chargino sector is identical 
to that of the MSSM. In the weak basis 



5, VTs, hi, h2, s 



the Lagrangian describing the neutralino mass matrix in the MNSSM is given by 

1 



^n-Mn^O + h.C. 



where 



Mr 





-MzSnjCp 






Mzc^Cfi 




-Mzs^Cfs 
Mzc^cp 




Vs 



-/i 



Mzs^sp 
-Mzc^sp 



V 



Vs 










(3.41) 



(3.42) 



(3.43) 
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with Cu, = ^l- sl = Mw/Mz. In Eq. ( glilD , B and VTs are the U(l)y and SU(2)i neutral 
gauginos, respectively, and hi, h2 and s are the corresponding Higgsino states of the chiral 
multiplets Hi, H2 and S. 

The neutralino mass matrix of the MNSSM given in Eq. ( |3.43| ) predicts a relatively 
light state, with mass smaller than 70 GeV. Since the neutralino mass matrix is identical to 
that of the PQ-symmetric extension of the MSSM, we call this light state axino a. The axino 
is predominantly a singlet field for values of the yU-parameter in the phenomenologically 
relevant range, i.e. for > 120 GeV. In order to have a first estimate of the axino mass, 
we assume that the gaugino mass parameters m~ and are very large, e.g. of order 500 
GeV and higher, such that the gauginos B and W3 decouple practically from the neutralino 
sector. The reduced 3x3 Higgsino-mass matrix, which is expressed in the subspace spanned 
by hi, h2 and s, can then be expanded in terms of v/vs-, thus yielding the axino mass 



V 



2 , 2A2 M2 



m-a ^ ^|/isin2/3| = — ^|sin2/3|. (3.44) 

This last formula proves to be a good approximation for \fj^\ ^ 200 GeV. 

There are strict collider limits on the axino-related parameters, which come from 
LEP2 and especially from the invisible width of the Z boson [|T^, in which case a new 
invisible decay channel for the Z boson into axino pairs opens up kinematically when < 
Mz/2. Assuming that the gauginos are decoupled from the neutralino mass matrix J^q, 
we find numerically that the axino mass is smaller than 45 GeV for values of ^ 150 GeV 
and A ~ (7^ ~ 0.65. Of course, such a numerical estimate crucially depends on the values 
of and For example, for relatively low values of mg and in the range 200-300 
GeV, gaugino-Higgsino mixing effects can no longer be neglected, and the upper limit on 
the /i-parameter is estimated to increase by 40-50 GeV. 

On the other hand, the upper bound on the branching ratio of the Z-boson invisible 
width due to a new-physics decay mode imposes the constraint 



B{Z^M) = ^^\g,,z\' < l.xlO-^ (3.45) 

at the 90% CL, where = g^j (47r) is the SU(2)/, weak fine structure constant and Vz = 
2.49 GeV is the total width of the Z boson. Moreover, in the seesaw-type approximation 
the aaZ -coupling is readily found to be 

9.az - J (4 -4) = ^#(4-4)- (3-46) 
The constraint in Eq. ( p.45| ), together with Eq. ( |3.46|) , leads to 

^^|cos2/5| < 0.122 . (3.47) 
9i ^^ 
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This last inequality can be translated into the following bound on the /i-parameter: 

\fx\ > 250 GeV, (3.48) 

for X ^ Qw and tan /9 ~ 2. The above exercise shows that in the MNSSM the LEP limits on 
the Z-boson invisible width give rise to a new exclusion range of fi values: 200 ^ J: 250 
GeV, for A ~ 0.65. However, this additional exclusion range of /i exhibits a quadratic 
dependence on A and completely disappears for values of A < 0.45. 



3.4 Effective Higgs-boson couplings 

Apart from the Higgs-boson masses, the effective couplings of the CP-even and CP-odd 
Higgs scalars to the and Z bosons are very essential for our phenomenological discussion 
in Section 5. These effective couplings are given by the interaction Lagrangians 

3 1 

£hvv = g^MwY.9H,vv{H^W;^W-'^' + —H^Z^Z^'), (3.49) 



Chaz = f^Y.llmA,z{H,d,A,)Z^, (3.50) 

^ — ^ < — 

where 9^ = 9^ — and 

QH^vv = cpO^, + spOg, (3.51) 

gu.A.z = Ot,{cpOi - s^OZ). (3.52) 

Here, we wish to remind the reader that the orthogonal matrix [O^) is related to the 
mixing of the CP-odd (CP-even) scalars and is defined in Eq. ( p.32|) ( (|3.40|) ). 



i=l '^w 
3 2 



It is now worth remarking that the effective couplings HiVV (with V = Z,W ) and 
HiAjZ satisfy the unitarity relations 



3 3 2 

9%vv = 1 > EE 9%A,z = 1 • (3-53) 

i=l 1=1 j=l 

In particular, in the limit in which A2 and H3 decouple as singlets, which is obtained 
for large \ts\ with the remaining parameters kept fixed, one recovers the known MSSM 



complementarity relations among the effective Higgs-boson couplings ||16|| : 



al^vv = al^A.z and gjj^yy = g^^A^z ■ (3-54) 

As an obvious consequence of the above decoupling limit, all couplings of the heavy Higgs 
scalars A2 and H3 to the W and Z bosons go to zero. 
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Another very important relation which involves the CP-even Higgs-boson masses and 
the respective couplings to the W and Z bosons is 



E gkvv Mm 



i=l 



4{Ml)n + 2sf3Cp{Ml)i2 + s}{Ml)22 
2A2 



Ml i cos2 2/3 + 



9i + 9 



12 



sin^ 2(3 
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167r2 



In 



2 2 



+ 



2X2 
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mi 



2 
t2 



In ^ 



.(3.55) 



This mass-coupling sum rule is very analogous to the one derived in pTI for the MSSM, 
where the RHS of Eq. ( p.55| ) is the squared lightest Higgs-boson mass in the decoupling 
limit of a heavy charged Higgs boson (see also Eq. (|3.17|) ). In this limit, only the Hi 
boson has non- vanishing couplings to the W and Z bosons |]T6[. As can be seen from 
Eq. ( p.55|) , the mass-coupling sum rule is independent of the charged Higgs boson mass, 
while it only weakly depends on fi at the one- loop order. The relations (|3.53|) and (|3.55|) , 
which are obviously valid for the case of the NMSSM as well, are very useful to reduce the 
number of independent effective Higgs-boson couplings and so achieve a better control on 
the numerical predictions for the Higgs-boson masses and couplings. 



4 MNSSM versus NMSSM 

Here, we shall compare the generic predictions for the Higgs-boson mass spectrum in the 
NMSSM, which includes the cubic singlet-superfield coupling, with those obtained in the 
MNSSM. For this purpose, we shall only focus on the tree-level structure of the Higgs sector 
of the NMSSM, as the dominant stop/top-radiative effects are identical for both models 
and have already been computed in Section 3.2. 

The often-discussed NMSSM is based on the Higgs superpotential 

Wffiggs = XSH^tT2H2 + (4.1) 

As usual, the complete tree-level Higgs potential is obtained by adding the relevant F- and 
D-teim contributions to the soft SUSY-breaking terms induced by the superpotential: 

- C^y = m2$l$i + m2$t$2 + mlS*S + [XA^S^\^2 + h.c.) + {^A^S^ + h.c.) 

-Ai($I$i)2 - A2($^$2)' - X3{<l>\<^i){H'^2) - (A4-A2)(<1'I<1'2)($5$i) 
+ S*S(<^\<^i + $^'^'2) + K^{S*Sf + [Ak5*2($|$2) + h.c.l . (4.2) 
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Furthermore, the minimization conditions are determined by requiring that the following 
tadpole parameters vanish: 
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(4.3) 



(4.4) 



— Vs rria + XA 
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K Vo 



(4.5) 



We should remark again that spontaneous CP violation is absent in the NMSSM at the 
tree level 0,^- Also, CP appears to be still a good symmetry of the NMSSM, even if 
(CP-conserving) large radiative stop-mixing effects were to be taken into account ||22|| . 

Considering the vanishing of the tadpole parameters given in Eqs. ( [4.3|) -( ^75D , it is 
not difficult to compute the charged Higgs-boson mass, and the CP-odd and CP-even mass 
matrices. More explicitly, the squared charged Higgs-boson mass is given by 
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2(0) 
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+ M, 
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(4.6) 



where the would-be /x-parameter is defined in Eq. ( |3.11|) . The entries of the tree-level 



CP-odd mass matrix Mp are found to be 
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3jfi j + 3-/uAk, 



(4.7) 



where Mm jg given by Eq. (|3.13| ). Finally, the entries of CP-even mass matrix Mg^^^ read 



(Mf))i2 = {Mf%i = -s,c,{Mm + M| - A^.;^) 

--(4c, Mm _ 2c,^,' - 



(Mf))3i 



Vs 



{Mf% = slMl + clMm^ 
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2 



{mT% = ^w(wM5°) + ^^^) - + 4^/i^ (4.8) 

From the above analytic expressions for Mfjl\ and the CP-odd and CP -even Higgs-boson 
mass matrices, Mp''-' and Mg^^\ it is now evident that the MSSM limit of the NMSSM is 
obtained for n, A 0, while holding k/A, fi, Ax and fixed. 

Parenthetically, we should remark that the Higgsino sector of the NMSSM is also 
different from the corresponding one in the MNSSM. Because of the presence of the operator 
I 5*^ in the superpotential ( [4.1D , the {55}-matrix element of the neutralino mass matrix in 
Eq. ( p.43|) receives the additional contribution: 

iMoh5 = -2^/^- (4.9) 

Note that if (A^o)55 < with /i < Oj^ this additional contribution to the predominantly 
singlet state in the NMSSM is constructive, rendering its mass larger than the axino mass 
in the MNSSM. However, for small positive values of k, e.g. k, < 0.1, and \^\ < 200 GeV, 
with yU < and A ~ 0.65, the {55}-matrix element {Aio)55 is positive and its contribution 
to the would-be axino mass is destructive, leading to light singlet masses smaller than rria- 

It is now important to notice that unlike the MNSSM case, the decoupled CP-even 
and CP-odd scalar singlets are no longer degenerate in the MSSM limit of the NMSSM. 
This fact is a manifestation of the violation of the mass sum rule ( p.l9|) in the case of the 
NMSSM. Specifically, in the NMSSM we find that 

i=l ^=l - ^4 A 

It is obvious that the mass sum rule ( |3.19| ) can be sizeably violated in the NMSSM for 
relatively large values of \k\ and or \Ai^\. In such cases, the violation of the mass sum 
rule becomes much larger than the one caused by radiative stop effects. 

The analytic expressions of the Higgs-boson masses in the NMSSM coincide with 

those of the MNSSM only in the PQ-symmetric limit, where k/X, ts, 0. Although 

this limit is unphysical as it leads to a theory with a visible axion, its vicinity, however, 

could define an acceptable region of parameter space where the predictions of the two 

models exhibit reasonable agreement. 

^Our choice of a negative /z-parameter is mainly dictated by the fact that 6 — > 57 imposes a stronger 
lower limit on positive values of fi for relatively small charged Higgs-boson masses, close to the present 
experimental bound, i.e. for Mh+ ^ 80 GeV Instead, for negative values of fi, the bound on fi can be 
dramatically relaxed up to the present LEP2 limit: > 90 GeV [Ea. 



E Mr - t Mlf' - MJ = 4 ^ ^= ( ^ + ^ - ^ ) . (4.10) 
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An interesting property of the tree-level CP-even mass matrix M^*-""* in the PQ- 
symmetric limit is that the interval of the allowed values is rather small. This interval 
may be determined by requiring that the determinant of Mg^^\ 
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with = m5°^ + Ml - IX^v"^ = mJ^? + M| - M^, be positive. Neglecting terms 
proportional to (|A^f^ — M|) cos^2/5 next to in Eq. (|4.11| ), we may approximate the 
determinant det(M^*-''-*) as 

det(M 



^2(0) ^ 
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Afi' - 2sin2 2/?M5°V' + ^ sin^ 2/3 M^W ( MfW 



Requiring now that det(Mg^°'') be positive gives the allowed fi^ interval: 
^sin22/3M2W (1-5) < /i^ < Isin22/3M5°) (1 + 5) 

with 
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Here, it is understood that 5 < 1 or, equivalently, M^+'' > M^. Especially for M^+'' = M^, 
for which the terms proportional to (|A^f^ — M|)cos^2/3 are no longer negligible with 
respect to = M|, the allowed range of /i^ becomes 
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Further insight into the predictions of the PQ-symmetric limit may be gained by analyzing 
the kinematic situation where M^+^ ^ (i.e. 5 -C 1) and /i^ = /i^jj = s^c^Mf*^°) which 
is approximately the middle point of the allowed /i^-interval. In this case, we obtain (see 
also Appendix B) 
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Thus, the hght Hi and Ai scalars decouple from the gauge bosons, whilst the H2 boson 
couples maximally to them with SM strength. Moreover, according to the mass-coupling 
sum rule ( |3.55| ), the if2-boson mass saturates the mass upper bound obeyed by the SM-like 
Higgs boson. Given that the length of the allowed /x^-interval is very small relative to yU^j^ 
for 5 <^ 1, one does not expect serious changes regarding the heaviest Higgs-boson masses 
Mhs, and Ma^ and the qualitative features of the Higgs to gauge-boson couplings as /i takes 
all other allowed values. 

A minimal deviation from the PQ-symmetric limit, in which the NMSSM could easily 
be compared with the MNSSM, is the limit k —>■ 0, with A, fi, Ax and held fixed. 
In fact, in this limit, the coupling A could be the largest, thereby allowing for the largest 
possible value for the lightest Higgs-boson mass Mh^- By the same token, the unwanted 
U(l)pQ symmetry gets broken by the trilinear soft SUSY-breaking self-coupling kA,^ of the 
singlet S. A corresponding limit of the MNSSM, which has the same number of independent 
parameters as in the NMSSM, is the one with ~^ 0; but A, /i. Ax and ts fixed. 
We should also bear in mind that vanishing of ts entails vanishing of mfg as well which 
eliminates the possibility of 0, with fixed. In this way, we compare essentially 

two models which only differ in soft operators of dimensionality d < 3. An additional 
reason that renders such a comparison very interesting is the fact that the dimensionful 
parameters, such as A^ and ts, remain unconstrained by perturbativity arguments, and 
hence could severely affect the structure of the mass matrices. 

The aforementioned physical limit allows for more direct comparisons of the NMSSM 
with the MNSSM. Equating the tadpole parameters T(^j, and T^^ pertinent to the two 
models yields, in this limit, the simple relation 

^ = (4.18) 
H A 

Moreover, in the same limit, except for (Mp°'*)22 and (Mg^^^)^^, all other elements of the 
mass matrices coincide as well. In the MNSSM Xts/n enters the elements (Mp '''*)22 and 
(M|*-°^)33 in exactly the same way, whereas in the NMSSM the corresponding parameter 
{K,fiAt^)/X appears in these two matrix elements with different coefficients, and most im- 
portantly, with different signs reflecting the violation of the tree- level mass sum rule (PTTPf). 



The fact that the determinants of the CP-odd Higgs-boson mass matrices in the MNSSM 
and NMSSM are proportional to Xts/fi and {K,fiAK)/X necessitates that the two parameters 
must be both positive. As a result of this, the matrix element (Mg^^^)^^ will be enhanced 
in the MNSSM, but reduced in the NMSSM. In addition, it is not difficult to see that the 
determinant of the tree-level CP-even mass matrix, det{Ms^^^), is a monotonically increas- 
ing function of {Mg^^^)s3 and is already negative if (M|*-'^^)33 vanishes. This last property 
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relies on the fact that the upper-left 2x2 submatrix is positive definite. Therefore, the 
larger the mass M^^ of the lightest CP-odd scalar is the larger (smaller) det(M^*-'^'*) is in the 
MNSSM (NMSSM). On the other hand, a very small value for Xts/fJ^ or {KfiAK)/X, which 
amounts to having a very light Ai, does not seriously affect Mg^^\ and hence no essential 
difference in the predictions for the Higgs spectrum between the MNSSM and NMSSM can 
be observed. In this region, both models are close to the PQ-symmetric limit. However, 
the difference between the two models becomes appreciable, once the parameters Xts/fi 
and {k,^Ai^)/X become large. The first parameter has no upper bound, whereas the second 
one is limited by the fact that {Mg^^^)s3 should be positive. Thus, only in the MNSSM 
case a significant departure from the PQ-symmetric limit is possible, which may change 
the situation drastically. For example, a //-independent contribution to (M|*-°^)33, say T^, 
changes the coefficient of /i^ in the expression ( PD for det(M|^°^), and as a consequence. 



the allowed interval for /i can now expand (or further shrink) for T positive (negative) 

f2( 



For the particular case that M'fj+ = M^, the interval of increases (decreases) by a 



factor 1 + T^/(A^t'^ sin^ 2/3). As we just observed, such an unconstrained (constrained) 
positive (negative) contribution is available in the MNSSM (NMSSM), i.e. = Xts/fi 
(T^ = —{nfiAf^)/X), where ts (Ak) should be regarded as a /i-dependent parameter. 

At this point, it should be stressed that our discussion of the NMSSM in the limit 
K —>■ 0, with A, fi, Ax and kA^ being kept fixed, by no means exhausts all possible 
predictions that the model offers for viable scenarios. Being close to the above limit requires 
that \k,/X\ <^ 1 <^ lA^Z/il. However, it is possible to considerably depart from this limit, 
even if |/t/A| is very small but non-zero, while avoiding the known problem associated with 
the presence of a visible axion. In order to better investigate alternative scenarios that 
avoid the presence of a visible axion, we compute the exact determinant of the CP-odd 
Higgs-scalar mass matrix 

det{Mf^) = 3^(/iA, - ^sin2/?AV - S^^Vj^^W. (4.19) 

We shall now examine other possible deviations from the PQ-symmetric limit, for which the 
fi values, however, are not very different from those determined by the allowed fi interval 
in Eq. ( |4.13| ). Under this assumption and the fact that \n/X\ is considered to be adequately 



small, the third term on the RHS of Eq. ( |4.19|) remains always subdominant; it actually 
diminishes the size of det(Mp°^) irrespectively of the sign of n/X. For ^ 1, the first 

term fiA^ on the RHS of Eq. ([4.19|) becomes dominant. In this case, this term should have 
the same sign as k/A, in compliance with our earlier requirement that (^/iA^) /A be positive. 
However, as I^k//^! is getting smaller, the second term |sin2/?A^f^ within the parentheses 
on the RHS of Eq. ( [4.19|) will then start playing an important role. For k/A>0(k/A<0), 
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this second term provides a lower (upper) bound on /iA^, which should not be saturated. 
In fact, the mass of the lightest CP-odd scalar depends crucially on the difference between 
these two first terms on the RHS of Eq. ( 4.19| ). It is then obvious that if n/\ is negative, 



yU^K could be negative, zero or even a positive quantity which is bounded from above. 

Having gained some insight from the above discussion, let us now consider the most 
general case without resorting to specific assumptions or kinematic approximations. Then, 
the requirement that det(Mp°'') in Eq. ( |4.19| ) be positive implies the constraint 

/i_ < < /i+ (4.20) 

with 



^l± = (6-5 



± sign(-) \/Al - 9sin2/352KAt;2 



(4.2i; 



where Af. > 9 sin 2/? S'^kXv'^. We see again that = is only allowed for k,/X < 0. In this 
case, /i^ is constrained to be in the range: 

< fi^ < -— sin2/3M5°^ (4.22) 

Here, it is also important to reiterate the fact that the requirement for a positive det (Mp"'*) 
constrains by itself the {33}-element of M|*^°'': 



{mT% < |w(.,c,m5°) - 2^/.^) + (4 - 36') ^/.^ (4.23) 



The constraint in Eq. ( 4.23 ) seems to favour negative values of k/A, as the upper limit on 
(M|*-°'')33 gets larger in this case. Furthermore, saturation of the upper bound in Eq. (^4.23|) 



leads to M^f^ = 0. 

As the key parameter |k/A| increases, the situation is getting more involved since 
new terms start playing a role. In particular, a term which deserves special attention is the 
one proportional to k^/i^/A^ that occurs in Eq. ( [i.23| ). This term becomes very important 



for larger values of Mjj+ which lead to smaller values of 6 and to larger values of jj'^ in 
accordance with Eqs. ( 4.13 ) and ( [4.14| ). In such a case, we may hope for an enlargement of 



the allowed interval of values, for which det(M^*-°'') is positive. Therefore, it would be 
interesting to investigate to which extent such a situation can indeed be realized, especially 
for low values of M^+^ for which 6 is not very small and {M'^^'^^)^^ appears to be more 



severely constrained. To this end, we shall consider the special case where M^^+^ = Mj 
i.e. 5 = 1. Then, after taking into account the constraint in Eq. ( |4.23| ) and making use 
of the fact that det{Mg^^^) increases monotonically with (M^^°^)33, the following inequality 
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may be derived: 
det(Ms 



^2(0). 
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(4.24) 

Assuming that the corresponding upper bound in Eq. (^^) is saturated (i.e. M^f^ = 0), 
then Eq. (|]2|) and the fact that det(M5^°^) > lead to 



< < - sin^ ' 
^ 2 



2/3 M^W 1 + - sin 2/3 + 
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cos^ 2/3 



(4.25) 



It is easy to see that for |k/A| <^ 1, the double inequahty in Eq. ([4.25|) reduces to our 
previous result found in Eq. ( |4.15| ). We observe now that for k/A > 0, the allowed interval 
of /z^ given by Eq. ( 4.25 ) shrinks as |k/A| increases. Instead, for k/A < with |k/A| 
increasing, the allowed interval gets larger and, especially for values of |k/A| close to unity, 
it may even become infinitely large. Of course, at this critical kinematic region, radiative 
corrections are expected to play the dominant role. Notwithstanding this fact, our tree- 
level results should still be indicative of the various tendencies which govern the kinematic 
parameters of the theory. As we will see below, however, values of |k/A| ~ 1 are not 
compatible with the largest possible value for A and hence with the largest value of the 
lightest Higgs-boson mass Mh^ . 

The Yukawa-type couplings n and A cannot be arbitrarily large, if we wish to preserve 
the good property of SUSY that perturbation theory be applicable up to the gauge unifi- 
cation scale Mu ~ 10^^ GeV [^. Therefore, upper limits on |A| and can be obtained 
by studying their renormalization-group (RG) evolution along with the corresponding ones 
of the strong coupling constant Qs and the t -quark Yukawa coupling ht |^|19],^ : 
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where t = \n{Q'^/M^). In writing the RG equations ( |4.2(j| ), we have ignored possible 
mass threshold effects of the SUSY particles while running from the t-quark-pole mass 
Mt = 175 GeV up to Mu ~ 10^^ GeV. In the RG analysis, we use the value for the strong 
fine-structure constant ^^(Mj) = g'^{Mt) / [Air) ^ 0.109. Furthermore, the running t-quark 
Yukawa coupling hf is determined by 

mt{Mt) 



ht{Mt) 



(4.27) 
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where v{Mt) = 174.1 GeV and 

Mt 

MMt) = . ^ 4 (4-28) 

is the known relation between the on-shell MS mass and Mt. For 3 ^ tan/5(Mi) < 10, we 
find the approximate upper bounds 

|A(Mi)| < 0.70, 0.63, 0.57, 0.44, 0.22, for \K{Mt)\ = 0, 0.3, 0.4, 0.5 and 0.6, (4.29) 

respectively. Correspondingly, for tan/5 ~ 2, we obtain 

\X{Mt)\ < 0.65, 0.59, 0.54, 0.42, 0.21. (4.30) 

The results in Eqs. ( ^.29| ) and ( ^.30| ) are in good agreement with Ref. [0J2^. From the 
above analysis, it is obvious that the largest value for |A(M()| is more naturally attained 
in the MNSSM (corresponding to K,{Mt) = 0) rather than in the NMSSM, as one would 
generically expect X{Mt) ~ K,{Mt) 7^ 0. This is another important difference between these 
two models. In particular, this implies that the MNSSM generically predicts higher masses 
for the lightest Higgs boson than the NMSSM. 

In the next section, we shall study the Higgs sector of the MNSSM more quantitatively 
and also compare our numerical predictions with those obtained in the NMSSM. 



5 Phenomenological discussion 

In this section, we shall discuss the phenomenology of the Higgs bosons in the MNSSM, 
and make comparisons of our predictions with those obtained in the NMSSM. 

At LEP2, the CP-even and CP-odd Higgs scalars, -f^i,2,3 and 2, are mainly produced 
through the Higgs-strahlung process e'^e'^ —>■ Z* —>■ ZHi or in pairs via e+e^ ^ Z* ^ HiAj. 
Analogous Higgs-boson production mechanisms can take place at Fermilab, where instead of 



electrons the initial states are the u and d quarks at the quark-parton level ^7j. Therefore, 
the necessary ingredients for our numerical discussion following below are the analytic 
expressions for the radiatively-corrected Higgs-boson masses and the effective Higgs-boson 
couplings to the gauge bosons. These analytic expressions pertaining to the MNSSM and 
NMSSM have been presented in Sections 3 and 4, respectively. 

There are several possible combinations in choosing the independent kinematic param- 
eters for the two supersymmetric extensions of the MSSM, the MNSSM and the NMSSM. 
For definiteness, for the MNSSM case, we consider 

tf3, M|f+ , /i. A, — and mj^, (5.1) 
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as free phenomenological parameters of the Higgs sector. As for the NMSSM, we take as 
input parameters 

tp , Mff+ , /i , A , K and ^ ^ . (5.2) 

A 

For both SUSY models, the stop-related parameters are chosen to have the typical values: 

Mq = Mt = 0.5 TeV, A = ITeV. (5.3) 

Here, we should remark that Eq. ( p.l| ) could in principle be absent, without spoiling 

the renormalizability of the theory. In this case, the U(1)pq symmetry of the MNSSM gets 
broken explicitly by the effectively generated tadpole parameter ts, which is a term of the 
lowest possible dimension, namely of dimension 1. Such a reduction of the renormalizable 
parameters is not possible in the NMSSM because of the presence of which violates 
U(l)pQ hardly. Therefore, it is fair to conclude that in the admissible limit mfj — ^ 0, the 
MNSSM under investigation represents the most economic, renormalizable scenario among 
the proposed non-minimal supersymmetric standard models. 

In Fig. ^ we display the dependence of the lightest Higgs boson Hi in the MNSSM 
with = on the would-be /i-parameter, for different values of the charged Higgs-boson 
mass, i.e. for Mh+ = 0.1, 0.3, 0.7 and 1 TeV. In Fig. 0(a), we choose the tadpole-parameter 
value Xts/n = 1 TeV^. As we are interested in maximal values for the lightest Higgs-boson 
mass Mhi which occur for low values of tan j3, i.e. for tan /5 = 2, we consider the largest 
allowed coupling A = 0.65, for which the MNSSM stays perturbative up to the gauge 
unification scale My ~ 10^^ GeV (see also discussion in Section 4). As can be seen from 
Fig. 1(a), the iJi-boson mass varies between 120 and 145 GeV depending on Mh+ for a 
wide range of fi values, which is significantly larger than the current experimental lower 
bound of 113.3 GeV on the SM-type Higgs boson. Furthermore, we observe an asymmetry 
of order 5 GeV in Mhi for large 300 GeV between positive and negative values of fi. 
This is because stop-radiative effects on M^i get enhanced for larger values of the stop- 
mixing parameter \Xt\ = \At — yu/ tan/5| which obviously result from large negative values 
of /i, provided jX^/max (Mq, M^)] < ^6 (cf. Eq. (HD). 

In Fig. H(b) we consider a smaller value for the tadpole parameter, i.e. Xts/f^ = 
0.04 TeV^. As in Fig. 0(a), we present numerical estimates of Mh^ as a function of /i, for 
the same discrete values of the charged Higgs-boson mass: Mh+ = 0.1, 0.3, 0.7 and 1 TeV. 
We find that the allowed range of fi becomes much smaller, but the maximum values of Mh^ 
are still very close to those obtained in Fig. |^(a). Most interestingly, we observe that the 
maxima of Mh^^ are located at almost the same n values found for the tadpole parameter of 
1 TeV^; the maxima are practically independent of the tadpole parameter, for all relevant 
values of \ts/ ^ = 0.01-1 TeV^. This feature that the allowed range of /i values shrinks as 
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Xts/fJ' gets smaller is in good agreement with our discussion in Section 4 concerning the 
CP-even Higgs-boson mass matrix in the PQ-symmetric limit. Specifically, for small values 
of Xts/fi, the allowed |/i|-ranges can be accurately determined by Eqs. ( |4.13| ) and ( ^4.15| ). 
In particular, the mean values of the allowed |/i|-ranges, which are approximately given by 
S13C13MH+ and are almost independent of Xts/fi', appear to describe well the location of the 
maxima of M^i- 

It is now very interesting to analyze a scenario within the context of the MNSSM, in 
which the charged Higgs boson has a relatively low mass, in the range Mh+ = 80-160 GeV, 
and may be accessed in next-round experiments at LEP2 and/or at the upgraded Tevatron 
collider. For this purpose, in Fig. |^ we display numerical estimates of the two lightest Higgs 
bosons Hi and H2 and their corresponding squared couplings to the Z boson as functions 
of the parameter fi, for Mfj+ = 80, 120 and 160 GeV. The other kinematic parameters 
are chosen to be the same as in Fig. H(a): tan/5 = 2, A = 0.65 and Xts/ ^ = 1 TeV^. Let 
us first consider the lowest experimentally allowed value for the charged Higgs-boson mass 
Mh+ = 80 GeV [|T],|rS|- Then, in Fig. |^(a) we notice that the i^i-boson mass cannot become 
larger than 105 GeV, whilst the next-to-lightest H2 boson can be as heavy as 146 GeV. 
As can be seen from Fig. ^(b), such a scenario is not excluded experimentally, since the 
HiZZ-coupling gets suppressed, i.e. gjj^zz ~ O-^, for |/^| ~ 350 GeV. In this scenario, the H2 
boson becomes SM type {H2 = Hsm), and is much heavier than H~^. This novel prediction 
of the MNSSM for viable scenarios with Mh+ < Mh^^ and low-values of tan/5 < 5 cannot 
be realized within the MSSM, even if CP-violating loop effects are included in the Higgs 
sector of the MSSM p8| , p!8[ .p] In fact, as we will see later on, neither the NMSSM can 
naturally accommodate scenarios with Mh+ ~ ^h^m ■> experimentally allowed values 

of > 90 GeV 0. 

As the charged Higgs boson becomes heavier in the MNSSM, the Hi boson also 
gets heavier and resembles the SM Higgs boson ifsM- Thus, from Fig. |^ we see that for 
Mh+ = 120 GeV, Mh, ^ 132 GeV, with g^^zz ~ 0.5, while for Mh+ = 160 GeV, it is 
Mh, < 142 GeV, with gj, 

jzz ~ 1- Furthermore, as we have already discussed in Section 
3.1, for the considered values of Mh+ much smaller than Xts/fi, the Higgs states A2 and 
H3 decouple and are almost degenerate with M^^ ^ Mjj,^ ^ Xts/fi- In particular, our 
numerical estimates confirm the relations: M^^ ^ ~ + 1-^^"^^ (^f- ■'^l- (P-16|) ), and 
9h^zz ~ 9h2A^z and gjj^zz ~ 9h^a^z (cf- Eq. ( |3.54|) ), which are only vahd in the above 
specific decoupling regime of the A2 and H^ bosons in the MNSSM. From Fig. ^(a), we see 
finally that for Mh+ = 160 GeV, the if2-boson mass is nearly /i-independent and equals 

^Using the code cph |29 based on one finds that only for extreme values of ^ 5 TeV and for 
tan/3 > 20, such a scenario might be made viable |]30| . 
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the Ai-boson mass Ma^ ~ 179 GeV. This result is just a consequence of the expected 
decouphng property of a heavy charged Higgs boson, with Mh+ = 160 GeV ^ Mz- 

In Fig. § we display predicted values for Mh^ and Mj^^, as well as for gfj^zz 
9h2ZZ the MNSSM, using the same input parameters as in Fig. but with tan/5 = 20, 
i.e. A = 0.65 and Xts/fi = 1 TeV^. We encounter a functional dependence of the evaluated 
kinematic parameters on fi, for Mh+ = 80, 120 and 160 GeV, which is qualitatively similar 
to the one presented in Fig. |^. Again, we see that the charged Higgs boson can be 
lighter than the Hi boson, even for large values of tan/5. Yet, we observe that for larger 
-ff^-boson masses, the squared ifi-boson coupling to the Z boson, g'jj^zz^ S^^s more rapidly 
to unity than in the tan (3 = 2 case. Here, we should emphasize again that as Mh+ becomes 
much larger than Mz, the predicted values for Mfji approach the one given by the square 
root of the RHS of Eq. ( p. 551) , where, of course, the term proportional to sin^ 2/3 is 
negligible. Therefore, only in this kinematic regime where both tan (3 and Mh+ are large, 
the predictions of the MNSSM will coincide with those of the MSSM. 

In our numerical analysis in connection with Fig. |^, we have already observed that 
for large values of Xts/fi but small values of Mh+, the Hi boson does not couple strongly 
to the Z boson, but it is rather the H2 boson which is SM-type. Actually, this kind 
of behaviour is encountered even for larger values of Mh+, provided Xts/fJ' is sufficiently 
small. As was already discussed in Section 4, the latter reflects the fact that the model 
approaches the PQ-symmetric limit in this case. In Fig. ^ we present numerically the 
//-dependence of the two lightest CP-even Higgs-boson masses, Mh^ and Mh2, and their 
respective couplings to the Z boson, gjj^zz ^"^^ Qh^zzi i'^ the MNSSM with = O5 
M}i^ = 200, 400, 600 and 800 GeV. In addition, we have selected the value of the tadpole 
parameter Xts/fJ' = 0.01328 TeV^. For this specific value of the tadpole parameter and 
for Mh+ = 400 GeV, we see that there is a value of /i where the Higgs states Hi and H2 
interchange their couplings to the Z boson, while being nearly degenerate having a mass 
close to the upper bound of Mh^- We shall denote by fi^, this specific value of /i at which 
a level crossing in the couplings of Hi and H2 occurs. Thus, for values of |/i| smaller than 
it is Qh^zz > Qh^zzi while this inequality of the squared couplings gets inverted for 
> |/i*|. If we now consider smaller values for = 200 GeV for the chosen value of 

Xts/ fJ', we observe from Fig. || that as \fi\ grows, gjj^zz starts higher than gjj^zz^ the 
crossing point of these two squared couplings is before Mh^ reaches its highest value. If we 
now take larger values for Mh+, e.g. Mh+ = 600, 800 GeV, we see that Qh^zz starts again 
higher, becomes almost unity with Mh^ close to its largest allowed value, according to the 
mass-coupling sum rule (|3.55| ), and intersects Qh^zz a smaller value of Mh^. The very 
special value of Xts/fi, for a given value of Mh+, for which Mh^ and become equal 
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at the highest possible value for Mh^ and Qh^zz ~ ~ 0-5 should be regarded as a 

critical point. Generically speaking, for values of \ts/ lower than the one corresponding 
to the critical point, the H2 boson couples predominantly to the Z boson. Instead, if Xts/fi 
is higher than its critical value, it is then the Hi boson that couples with SM strength. In 
addition, in Fig. ^we see that almost independently of Mh+, the squared couplings Qh^zz 
and g\^zz remain comparable for a wide range of values. The latter is an indication of 
the fact that the critical value of At^/yU depends only weakly on the charged Higgs-boson 
mass Mu^ and has a value close to 0.01 TeV^, for 0.3 TeV Mn^^ < 1 TeV, where the 
remaining independent kinematic parameters are held fixed. 

We shall now analyze the predictions of the MNSSM for relatively small values of 
the tadpole parameter Xtsll^- In Fig. |], we display numerical estimates of M^^ and M^^, 
as well as of g\^zz Sh-^zz-- functions of the /i-parameter, for \tsl\i = 0.0026 TeV^. 
As for charged Higgs-boson masses, we choose Mh+ = 0.3, 0.5, 0.7 and 1 TeV. It is easy 
to see that, to a good approximation, the functional dependence of the masses of the two 
lightest CP-even Higgs bosons Hi and H2 are insensitive to the value of Mh+- In this 
scenario of the MNSSM, the H2 boson has always the strongest coupling to the Z boson. 
Although not displayed in Fig. the mass of the lightest CP-odd scalar Ai is found to be 
Ma^ ~ 50 GeV and is almost independent of Mh+. In addition, the CP-odd Higgs scalar 
Ai has suppressed couplings to the Z and Hi bosons, i.e. Qh^a^z ~ 10"^, and therefore 
can escape detection at LEP2. In Fig. |^ we notice finally that the allowed intervals of 
/i values become even shorter than those found in the previous scenarios of the MNSSM. 
These results are all consequences of our choice of a relatively small value for the tadpole 
parameter and are in good qualitative agreement with our discussion in Section 4 pertaining 
to the PQ-symmetric limit. 

It is very interesting to examine the consequences of the presence of a non- vanishing 
effective Fj-tadpole term on the Higgs-boson mass spectrum of the MNSSM. Therefore, 
in Fig. we plot the dependence of the CP-even Higgs-boson masses Mh^ and and 
the squared couplings g'jj^zz and g\i^zz-: as functions of the /^-parameter, for ts = —1 TeV^ 
and 771^2 = 0.325 TeV^. Because of the close relationship between ts and m^2) we are now 
compelled to treat ts as a yU-independent constant. In fact, for Msusy = 1 TeV, we can 
easily compute from Eq. ( |3^ ) that the adopted values for ts and correspond to the 
typical values of ^s and ^p: C,s = —1 and = 1/2- To enable a direct comparison with 
Fig. 1^, we choose the same values as in Fig. ^ for the remaining kinematic parameters of 
the theory. From Fig. |^, we see that the presence of a non-vanishing, positive tadpole 
term can shift the maxima of Mh^ and towards larger values of whereas 
all other features found in Fig. |^ are retained. In Fig. 0(a), we have also displayed the 
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dependence of the mass Ma^ of the hghtest CP-odd scalar Ai, as a function of /i. We see 
that decreases with decreasing. This kinematic behaviour originates from the fact 
that the contribution of the off-diagonal terms to the CP-odd mass matrix becomes rather 
significant for smaller values of \fi\. Instead, for larger values of the corresponding 
contribution of the off-diagonal terms is smaller, and leads to the mass relation M^^ ~ Ma- 

Unlike the MSSM, the charged Higgs-boson cannot be arbitrarily heavy in the 
MNSSM for fixed given values of tan /3 and A, and for natural choices of Xts/fJ' and the soft 
squark masses, i.e. for Xts/ fi, Mq, < 1 TeV^. Figure | displays the dependence of the 
maximum of the lightest Higgs-boson mass, max (Mf/J, as a function of Mh+, for tan (3 = 2, 
A = 0.65 and for two different values of the tadpole parameter: Xts/fi = 0.04 and 1 TeV^. 
The coupling of the Hi scalar to the Z boson becomes SM-type, for Mh+ <; 150 GeV. From 
Fig. H, it is then easy to see that the current LEP2 lower bound on max {Mh^ ) implies the 
approximate upper limit on Mh+- Mh+ < 2.7 TeV, almost independently of Xts/fi- This 
result may be understood as follows. As the mass of the charged Higgs boson increases, the 
maximum of occurs for larger values of which is a consequence of the tree-level 
structure of the CP-even Higgs-boson mass matrix in Eq. ( p.l5| ). On the other hand, the 
larger the value of becomes the larger the stop-mixing parameter \Xt\ = \At — fi/tpl 
is getting. Thus, when |X( |/max (Mq, M^) > y^, the contributions of the stop-radiative 
effects to the lightest Higgs-boson mass Mh^ become negative, and so drive max (Mf^J to 
unphysical values. For the very same reasons, a similar dependence of max [Mh^ ) on Mh+ 
is found to apply to the NMSSM case as well. 

For comparison, we shall now investigate a few representative scenarios within the 
context of the NMSSM. As a first example, we consider the scenario with tan (3 = 2, X = 
0.65, K = 0.01 and {KfiAf^)/X = 0.0026 TeV^. Figure ^exhibits the numerical predictions 
for the two lightest Higgs-boson masses, M^-^ and Mh2, and their corresponding squared 
couplings to the Z boson, gjj^zz ^"^^ 9hiZZ-- functions of the /^-parameter. We also vary 
the charged Higgs-boson mass in a discrete manner, i.e. Mii+ = 0.3, 0.5, 0.7 and 1 TeV. 
We observe that Mh^ and are practically independent of Mh+, with Mh^ consistently 
below 80 GeV. Such low values of Mh-^ are still acceptable at LEP2, in the range of /x 
values where g^n-^zz ~ 0.07. In this scenario, the H2 boson has a SM-type coupling to the 
Z boson. Also, the mass of the lighest CP-odd Higgs scalar Ma^ is almost independent of 
Mh+ and comes out to be slightly higher than Mh^. The NMSSM under discussion, with 
the chosen low value of k ^ 0.01, may be considered to adequately describe the limiting 
scenario where k — > and nA^^ is held fixed. This last fact enables one to directly compare 
the present scenario of the NMSSM with the MNSSM where the tadpole parameter Xts/fi 
is set to the same value with that of {KfiAi^)/X, i.e. Xts/fJ' = 0.0026 TeV^. Such a scenario 
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in connection with the MNSSM has already been analyzed above in Fig. ^. Thus, if we 
now compare Fig. ^ with Fig. ^, we observe resembling numerical predictions for the Higgs- 
boson masses and couplings in the two models. The only visible difference between them is 
that in the MNSSM, the lightest CP-even Higgs boson Hi is consistently 30 GeV heavier 
than the corresponding one in the NMSSM, while the mass M^^ of the lightest CP-odd 
scalar is about 30 GeV lower. These findings are in excellent agreement with our discussion 
in Section 4. 



We shall now analyze in Fig. |10| a second scenario of the NMSSM, in which the 
Yukawa-type coupling k is larger, but with the expression {K,fiAK,)/X being held fixed again, 
i.e. K = 0.1 and {K,fiAK,)/X = 0.0026 TeV^. In Fig. IC, we also vary the charged Higgs-boson 
mass in the same way as in Fig. |^: Mh+ = 0.3, 0.5, 0.7 and 1 TeV. In this scenario, the 
ratio |v4K//i| varies from 1.69 for = 100 GeV up to 0.096 for = 420 GeV, namely the 
ratio IAk//^! is no longer much larger than 1 for all relevant values of Furthermore, as 
|/i| increases, the strong inequality lAsZ/il ^ \i^/M S^ts gradually violated as well. As a 
consequence, as the charged Higgs-boson mass Mh+ takes higher values, the picture starts 
changing in comparison with Fig. y. To be precise, as Mh+ becomes larger, we observe a 
progressive enhancement of the maximum of the ifi-boson mass Mh^^ and of its respective 
squared coupling to the Z boson g'jj^zz'^ ^^e values of Mh^ and g\i^zz approach those of 
Mji^ and Qh^zz^ respectively. In particular, when Mh+ approaches 1 TeV, a level crossing 
effect in the masses and couplings of the Hi and H2 bosons takes place and the Hi boson 
becomes SM-type. In addition, the mass of the lightest CP-odd Higgs scalar Ma^ gets very 
small, i.e. Ma^ ~ 15 GeV, resulting from a partial cancellation of the first two terms on 
the RHS of Eq. (^191) . It is obvious that with increasing \k\ and the predictions of the 
NMSSM start slowly resembling those of the MNSSM with Xts/n being in the vicinity of 
its critical value. 



On the other hand, as the charged Higgs-boson mass decreases, we notice in Fig. 10 



that viable scenarios occur for smaller values of In fact, within the specific NMSSM 
with K, = 0.1 that we have been considering here, the experimental constraint, ^ 90 GeV 
T5| , implies that Mh+ cannot be lighter than 180 GeV. Of course, such a scenario could 



be directly excluded from the fact that for small positive values of k ~ 0.1, the lightest 
singlino state contributes significantly to the Z-boson invisible width (see also discussion 
after Eq. ( [4.9|) ). For this reason, we present in Fig. numerical estimates for a related 
scenario with negative k, i.e. n = —0.1. We also choose a smaller value for {KfiAi^)/X, i.e. 
(K/iA^)/A = -0.0021 TeV^ so as to obtain a light CP-odd Higgs state Ai. In Fig. |Tl|, we 
observe again the same characteristics as in Fig. 10, namely as the charged Higgs-boson 
mass decreases, viable scenarios take place for smaller values of leading to a similar 
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lower bound of about 180 GeV on Mh+- In fact, after having carefully explored all the 
relevant parameter space of the NMSSM with tan/? = 2, we found that this is a general 
feature of the NMSSM for any perturbative value of A and k (see also discussion below). 
Consequently, as in the MSSM, the SM-type Higgs boson in the NMSSM is also predicted 
to be lighter than the charged Higgs boson. 

To get a better understanding of this last phenomenological feature of the NMSSM, 
it is very instructive to analyze a scenario where the Yukawa-type couplings k, and A are 
comparable in size. Specifically, we choose A = 0.5 and k, = —0.45. According to our 
discussion in Section 4, the parameters of this model have been chosen in a way such that 
the charged Higgs boson might be allowed to become lighter than the one predicted in 
the previous scenarios of the NMSSM. Furthermore, in order to obtain the largest possible 
values for the masses of the CP-even Higgs scalars, we always fix by the requirement that 
the Ai boson be extremely light of the order of a few GeV.|U] Having the above in mind, we 
present in Fig. |l^ numerical estimates of Mh^ and Mh2, and g'jj^zz Qh^zzi functions 
of the /i-parameter, for charged Higgs-boson masses M/^+ = 120, 400 and 800 GeV. We 
observe that for |/i| > 100 GeV, the Hx boson is always SM-type. In addition, for M^i^ = 
120 GeV, the mass of the Hi boson has a maximum of ~ 113 GeV at = 100 GeV 
with Qh^zz — 0-5; which is close to the present experimental lower bound of LEP2 [|ll. In 
this scenario, the next-to-lightest CP-even Higgs boson H2 has a smaller coupling to the Z 
boson and its mass varies between 120-130 GeV. For larger values of Mh+, the Hi boson 
is always SM-type, with M^^ ~ 120-130 GeV for a wide range of values, whilst the H2 
boson is very heavy and decoupled from the lightest Higgs sector. 

Our numerical analysis as presented above in Fig. ^ explicitly demonstrates that for 
Mh+ = 120 GeV, the mass Mh^ of the lightest CP-even Higgs boson becomes acceptable 
only within a very narrow interval of \fi\, which is, however, close to its current lowest 
bound as set by LEP2 |15]. Thus, even within this optimized scenario of the NMSSM with 
|k/A| ~ 1, the Hi boson cannot become heavier than the charged Higgs boson. Therefore, 
we reach the conclusion that a possible discovery of a charged Higgs boson lighter than 
120-130 GeV and a SM-type Higgs boson heavier than 130-140 GeV can only be naturally 
accounted for within the MNSSM. 



"Despite its similarity, our scenario differs from the one discussed in |^ very recently. In our case, the 
tree-level values of A\ and required for Ma^ ~ 0, are not forced to be suppressed. The latter turns 
out to be the case only within a very narrow range of /x values close to the upper-end of the interval given 
by Eq. ( p2|) . 
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6 Conclusions 



We have considered the simplest extension of the minimal supersymmetric standard model, 
in which the /i-parameter has been promoted to a dynamical variable by means of a gauge- 
singlet superfield 5", with the linear, quadratic and cubic singlet-superfield terms, S, S"^ 
and S^, absent from the superpotential. Moreover, we have assumed that the breaking of 
SUSY in the observable sector is communicated gravitationally by a set of hidden-sector 
superfields which break = 1 supergravity spontaneously. In such a supergravity scenario, 
the absence of harmful destabilizing tadpole divergences at lower loop levels can be assured 
by forcing the complete superpotential and Kahler potential to respect specific discrete 
R symmetries. In particular, we have been able to show that with the imposition of the 
discrete R symmetries and Z^, the potentially dangerous tadpole divergences first 
appear at the six- and seven-loop levels, respectively, and hence are naturally suppressed 
to the order of the electroweak scale, without destabilizing the gauge hierarchy. 

The MNSSM we have been studying in this paper has a number of appealing field- 
theoretic and phenomenological features, which may be summarized as follows: 

• The model provides a natural solution to the so-called /x-problem of the MSSM, since 
the value of the /i-parameter can now be directly set by the VEV of the gauge-singlet 
superfield 5" which is of the required order of Msusy, as a consequence of the Z^ and 
Z^ symmetries. 

• The presence of the effectively generated tadpole terms linear in S and Fs (or S) 
breaks explicitly the continuous U(1)pq and its discrete subgroup Z^. Thus, the 
model offers a natural solution to the visible axion and cosmological domain-wall 
problems. 

• Depending on the underlying mechanism of SUSY breaking, the effective tadpole 
proportional to Fs could in principle be absent from the model. Such a reduction 
of the renormalizable operators does not thwart the renormalizability of the the- 
ory. The resulting renormalizable low-energy scenario has one parameter less than 
the frequently-discussed NMSSM with the cubic singlet-superfield term ^S^ present; 
it therefore represents the most economic, renormalizable version among the non- 
minimal supersymmetric models proposed in the literature. 



As opposed to the NMSSM, the MNSSM satisfies the tree- level mass sum rule ( ^.19 ) 



which is very analogous to the corresponding one of the MSSM [16|. This striking 



analogy to the MSSM allows us to advocate that the Higgs sector of the MNSSM 
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differs indeed minimally from the one of the MSSM, i.e. the introduced model truly 
constitutes the minimal supersymmetric extension of the MSSM. In the NMSSM, the 
violation of the mass sum rule can become much larger than the one induced by the 
one- loop stop/top effects, especially for relatively large values of \fi\ and \Ak\. 

A generic prediction of the non-minimal supersymmetric standard models is that for 
low values of tan (3, the lightest CP-even Higgs-boson mass M^^ increases significantly 
with growing |A| [cf. Eq. (3.18)]. Since in the MNSSM A can take its maximum 



allowed value naturally corresponding to the NMSSM with n = 0, the value of Mh^ is 
predicted to be the highest, after the dominant stop-loop effects have been included, 
i.e. Mh^ < 145 GeV. Therefore, such a scenario can only be decisively tested by the 
upgraded Tevatron collider at Fermilab and by the Large Hadron Collider (LHC) at 
CERN. 

The MNSSM can comfortably predict viable scenarios, where the mass of the charged 
Higgs boson H'^ is in the range: 80 GeV < Mfj+ < 3 TeV, for phenomenologically 
relevant values of ^ 90 GeV. In particular, numerical estimates in Section 5 reveal 
that a possible discovery of a charged Higgs boson, with Mh+ < 120 GeV, and a 
neutral Higgs boson, with M^i ~ 130 GeV, can only be naturally accounted for within 
the MNSSM, whereas the NMSSM would be highly disfavoured. This important 
phenomenological feature of the MNSSM, which is very helpful to discriminate it from 
the NMSSM, is a reflection of a new non-trivial decoupling limit due to a large tadpole 
\ts\, which is only attainable in the MNSSM (see also discussion of the paragraph 
that includes Eq. (|3.16| )). 



For scenarios with Mh+ ~ 200 GeV, the distinction between the MNSSM and the 
NMSSM becomes more difficult. In this case, additional experimental information 
would be necessary to distinguish the two SUSY extensions of the MSSM, resulting 
from a precise determination of the masses, the widths, the branching ratios and the 
production cross sections of the CP-even and CP-odd Higgs bosons. Nevertheless, 
if the tadpole parameter Xts/fJ' becomes much larger than with the remaining 
kinematic parameters held fixed, the Higgs states H3 and A2 will be predominantly 
singlets. As an important phenomenological consequence of this, the complemen- 
tarity relations ( p.54| ) between the H12ZZ- and H2^iAiZ- couplings will then hold 



approximately true in the MNSSM. However, these relations will be generically vio- 
lated in the NMSSM, as there is no analogous decoupling limit in the latter model, 
in which the states and A2 could decouple as singlets. 

The MNSSM also predicts the existence of a light neutralino, the axino. The axino 
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is predominantly a singlet field, for > 120 GeV. LEP limits on the Z-boson invisible 
width lead to the additional constraint: 200 < < 250 GeV, for A ~ 0.65. However, such 
a constraint disappears completely for smaller values of A, namely for A ^ 0.45. In fact, the 
axino may become the LSP in the MNSSM. In this paper we shall not address the issues 
associated with the cosmological consequences of the axino on the reheating temperature 



of the Universe and on the dark-matter problem. A detailed discussion of these issues 
may be given elsewhere. 

The present study has shown that the MNSSM is a viable scenario, which departs 
minimally from the MSSM, having a large number of appealing field-theoretic and phen- 
emenological features. Even though further refinements of our treatment of loop effects 
might be very useful, such as the inclusion of one-loop D-teim contributions to the effective 
potential and the computation of two-loop leading logarithmic corrections, our predictions 
for the Higgs-boson mass spectrum as well as the results of our comparative analysis be- 
tween the MNSSM studied here and the frequently-discussed NMSSM are not expected to 
modify dramatically. In particular, we find that the MNSSM can naturally predict viable 
scenarios in which the charged Higgs boson H'^ is much lighter than the neutral Higgs boson 
with SM-type coupling to the Z boson. The planned colliders, i.e. the upgraded Tevatron 
collider and the LHC, have the potential capabilities to test such interesting scenarios with 
a relatively light H^, as well as probe large domains of the Higgs-sector structure of this 
truly minimal supersymmetric extension of the MSSM, the MNSSM. 
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A Non-destabilizing tadpole divergences 



Employing standard power counting rules we shall show the absence of harmful 

tadpole divergences up to a sufficiently high loop order n, i.e. n < 5, within the context of 
the supergravity scenarios described in Section 2. 

It is useful to briefly review flrst the sufficient conditions that govern the absence of 
harmful tadpole divergences. To this end, let us consider a supergraph with one external 
leg, i.e. a tadpole graph. The tadpole graph may involve a number Vd of superpotential 
vertices of dimension d + 3, which are of the form z'^'^^/Mp where z represents a generic 
chiral superfleld, and a number Ud of Kahler-potential vertices of dimension d + 2, which 
have the form z'^^'^/Mp. Then, the superflcial degree of divergence of the tadpole graph, 
e.g. that of S, is given by 

D = 1 + Y.dVd + Y.dUd, (A.l) 

d d 

which leads to a contribution to the effective potential 

1 ^^'D+j:,dv,+j:,du, 

(A.2) 

where n counts the number of loops and Msusy is the soft SUSY-breaking scale. In ob- 
taining the last step of Eq. ( |A.2| ), we have used A ~ Mp as a natural energy cut-off scale. 
This very last step in Eq. ( [A. 21 ) shows that a tadpole contribution to the effective potential 
is proportional to one power of Mp at most. Such tadpole contributions which remain pro- 
portional to Mp will be referred to as "harmful" to be distinguished from the "harmless" 
ones in which the cut-off dependence disappears. In this context, an additional requirement 
for a tadpole graph to be harmful is that D be an even number. Finally, the degree of 
divergence can also be determined by the number of loops n and superpotential vertices Vd 
through the relation 

D = 2n -J2Vd. (A.3) 

d 

In summary, one flnds that a set of vertices produces a harmful tadpole divergence if the 
following equalities are all simultaneously satisfied: 

D = 1 + Y^dVd + Y.dUd = 2n - Y^Vd = even, (A.4) 

d d d 

with D>2. 

In the next two subsections, we shall apply the power counting rule of superficial 
divergences, stated in Eq. ( |A.4|) , to the two models based on the discrete i?-symmetries 
and Z^. 
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A.l The case 



Here, we shall show that the potentially harmful tadpole divergences are absent up to five 
loops. Alternatively, we shall prove that it is impossible to construct a tadpole diagram 



D < 10. Suppose now that at least one superpotential vertex is involved in a tadpole 
supergraph. Based on Eq. (|A.3|) , we see that we need at least two superpotential vertices 
to form a tadpole graph with D even, i.e. Yl,d^d > 2. Thus, for n = 5, one has D < 8, 
and by virtue of Eq. ( [A.l|) , it is Y.d.dVd < 7 and d <7. In the case that no superpotential 
vertices are involved, we have the relation D = 1 + J2d dUd < 10 or J2d dUd < 9 on account 
of Eq. ( [A. 11 ). We consider it obvious that it is impossible to form a tadpole graph with only 
one Kahler-potential vertex of c? = 9. This observation excludes Kahler-potential operators 
of c? = 9. Furthermore, as we will see below, the imposition of on the complete Kahler 
potential does not permit operators oid = 1. If we now wish to satisfy the above constraint 
Z^d dUd < 9 with two vertices, we then need one operator oid = 2 and another one of (i = 7; 
the latter is the Kahler-potential term of the highest dimensionality that could participate 
into a harmful divergent tadpole graph with n < 5. Consequently, we reach the conclusion 
that only superpotential and Kahler-potential vertices with d < 7 will be of relevance here. 

We shall confine ourselves to a minimal model, in which only the superfields Hi, H2 
and S are present and ignore quark and lepton superfields, as they do not couple directly 
to S; the inclusion of the fermion superfields is straightforward and does not alter our 
results. Moreover, we shall not include in the list of Kahler-potential terms those obtained 
by multiplying the latter with any power of hIHi, hIH2, S*S. The reason is that the 
omitted terms generate graphs of higher loop order than the included ones. 

Having the above in mind, we are now able to list all superpotential and Kahler- 
potential terms of d <7, respecting the discrete _R-symmetry (cf. Eq. ([^.6[)): 



from the sets of vertices which satisfy the condition ( |A.4D for n < 5, corresponding to 



W : Wf 




Wi = 




6(9) + h.c. 



(A.5) 



K^^^ = S*S, 
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K, = + h.c., K, = + h.c., 

Kr = j^, + h.c., (A.6) 

where 6{6) is the usual Grassmann- valued 5-function. Notice that the terms Kq^\ K^^ 

(3) 

and Kq represent the usual Higgs-superfield propagators and have no direct effect on 
our power counting rules. These terms are merely needed to contract the superfields in 
propagator lines and so form a loop supergraph. Furthermore, from Eqs. ( |A.5| ) and ( |A.6D , 
we observe that forbids the appearance of superpotential operators of = 2 {W2) and 
of Kahler-potential terms of c? = 1,4 {Ki^K^). 

In the following, we shall systematically analyze all possible sets of vertices compatible 
with the conditions in Eq. ( |A.4| ) up to five loops. At the one-loop level {n = 1), with 
J2d^d = 0, we readily find from Eq. ( |A. 1| ) that J^d^Ud = 1, entailing the absence of 
contributing operators. The situation becomes increasingly more involved for n = 2, 3, 4 
and 5. More explicitly, our systematic search for the existence of possible harmful tadpoles 
may be summarized as follows: 



I. n = 2 : 



II. n 



a) D = 2, = 2, Y.dVd + Y.dUd = I: 

d d d 

{Wo, Wi}; 

b) D = 4, = 0, Y.dUd = 3: 

{K?}- (A.7) 

a) D = 2, = 4, + T.dUd = I: 

d d d 

b) D = 4, E^rf = 2, Y.d^d + Y.dUd = ?>: 

d d d 

{Wo, W^}, {Wo, W^, K^}, {2Wo, Kf}- 

c) D = 6, Y.^d = 0, Y.dUd = 5: 

{K,}, {K2, (A.8) 
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D ^ 2, Y.Vd = 6, J2dVd + Y.dUd = 1 : 

d d d 

{5^0, Wi}; 

L> = 4, E^ci = 4, + ^dC/rf = 3: 

{3W^o, W^s}, Vo, 3W^i}, {3V0, W^i, ^2}, 

{41^0, Ki^y, 

D ^ 6, E^ci = 2, J^dV^, + J^dC/rf = 5: 

W, W5}, {Wi, 1^4}, W^3, {Wo, Wi, 2K2}, 

{2W,, 4'^}, {2Wo, K,}, {2Wo, K2, K^^h 
D = 8, Y^^d = 0, Y.dUd = 7: 

d d 

{2X2, }, {i^2, i^s}, {i^r}. (A.9) 
D ^ 2, E^ci = 8, ^dKi + ^dC/d = 1: 

d d d 

D ^ 4, = 6, Y.dVd + Y.dUd ^ 

d d d 

D = 6, = 4, E'^^ci + E'^^'i = 5: 

{31^0, W^s}, Wo, W^i, ^^4}!" W, 2M^i, W^a}, 
{31^0, W^3, K^}, W, 31^1, K2}, {21^0, SW^i, Xf}, 
{31^0, W^i, 2i^2}, {41^0, i^5}, {4V^o, i^2, i^?^}; 
D = 8, E^d = 2, E'^^ci + E'^^'i = 7: 

{W^l, W^4, i^2}, W, W^4, i^?^}, W, W^3, i^P}, 

{1^0, W^3, 27^2}, {21^1, i^5}, {2W^i, i^2, i^?}, 

{2Wo, Kr}, {2Wo, K^, K,}, {2Wo, 2i^2, 
D = 10, E^ci = 0, E^t^<i = 9: 

d d 
{i^2, i^7}, {2i^2, K,}, {3K2, 4'^}, {i^6, 4^}, 
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{kI\ K!f\ irf }. (A.IO) 

Here i,j,k = 1,2,3. The remaining task is to show that the sets of vertices hsted above 
do not produce tadpole graphs. This can be best verified case by case algebraically in 
the following manner. First, we multiply all the vertices belonging to a set and formally 
substitute H\Hi, H2H2 and S*S with 1 into the product of vertices. Then, we examine 
whether terms linear in S or S* survive in the resulting expression. In this way, we have 
carefully checked that there are no such terms linear in S or S* for all sets of vertices listed 
in Eqs. ( |A.7|) -( [A-.10 ), thus implying the absence of harmful tadpole graphs up to five loops. 



At a higher loop level, we can construct tadpole supergraphs by making free use of 
the renormalizable superpotential vertex Wq in Eq. ( [A. 51 ) together with some of the above 
vertices e.g., the higher-dimensional Kahler-potential vertices K2 and defined in Eq. 



( |A.6| ). Specifically, we find that the set of vertices 

{4Wo, K2, K,} (A.ll) 

leads to the typical six- loop tadpole graph depicted in Fig. |l](a). Also, it is not difficult 
to see that the above graph is actually a harmful divergent one since the set of vertices in 
( |A.11| ) satisfies the global constraint of Eq. (|A.4|) , with n = 6, D = 8 = even, J2d^d = 4, 



and Ed dVd + Ed dUd = Y.ddUd = 7. 



A. 2 The case 

In this section we shall show that the symmetry prohibits the presence of all possible 
harmful tadpole divergences up to six loops. Following a line of arguments similar to the 
case, we conclude that only superpotential and Kahler-potential operators with d < 9 
are of interest in this case. Therefore, we list all possible vertices of d < 9, respecting the 
discrete i?-symmetry Z^^ (cf. Eq. (|2.9|)): 

W : Wo = S{HliT2H2)5{e) + h.c, W2 = jpSie) + h.c, 

^ ^^5{e) + h.c. (A.12) 
K: K'i'^ = H\%, irf = #|#2, K'i^ = S*S, 
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- Ml. + ^' - Ml + 

Note that the symmetry forbids the occurrence of superpotential operators of c? = 1, 4, 6 
{Wi, Wi, Wq) as well as of Kahler-potential terms of ci = 1, 2, 7 (i^i, Kj). 

As we did for the case, we shall determine all possible sets of vertices compatible 
with the conditions in Eq. ( |A.4|) up to six loops. Again, it is not difficult to see that at the 
one-loop level (n = 1), with J2d^d = and J^d^Ud = 1, one is unable to find contributing 
operators. Furthermore, the absence of = 1 operators leads to the constraint 

Y^dVd + Y.dUd = D - I > 1, (A.14) 

d d 

i.e. D > 2. This observation simplifies further the search for the existence of possible 
harmful tadpoles. Thus, for n = 2, 3, 4, 5 and 6, we find the set of vertices: 



I. n = 2 : 



II. n = 3 : 



III. n = 4 : 



a) D = 4, ^V^rf = 0, J2dUd = 3: 

(A.15) 

a) D = 4, J^V^rf = 2, + Y.dUd = 

d d d 

{Wo, Ws}, {2Wo, Ki% 

b) D = Q, = 0' J2dUd = 5: 

{K,}. ' ' (A.16) 



a) D = A, = 4, J2dVd + Y.dUd = 3 

d d d 

{3Wo, W3}, {AWo, Ki% 

b) D = Q, J^^d = 2, J2dVd + Y.dUd = 5 

d d d 
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{Wo, Wr,}, {W2, W3}, {Wo, W2, }, {2Wo, K,}; 
D ^ 8, = 0, Y.dUd = 7: 

{Ki\ K,}. (A.17) 



= 4, E^ci = 6, Edyd + E^^d = 3: 

d d 

[hWo, VTa}, {6H^o, }; 

L> = 6, E^ci = 4, Y.dVd + Y.dUa ^ 

d d d 

D = 8, E^ci = 2, E^^^rf + E^^f^ci = 7: 

W, W^r}, {W2, W^}, {Wo, W2, K^}, 

{Wo, W,, K,}, {2W2, i^?}, {2Wo, kI\ K,}; 

D = 10, - 0, Y^dUd = 9: 



D ^ A, E^c^ = 8' Y.dVd + Y.dUd ^ 

d d d 

{IWo, W^}, {mo, }; 

D = = 6, Y.d^d + Y.dUd = 5 : 

d d d 

{^0, W^}, {AWo, W2, W^}, {^0, W2, }, {oWo, K^y, 

D = 8, E^rf = 4, E^^rf + T.dUd = 7: 
{3W^o, Wr}, \2Wo, W2, W^}', {Wo, 2W2, W^}, 
{2Wo, 2W2, i^?^}, {?>Wo, W^, {?>Wo, W2, K^}, 

{AWo, Kf, X4}; 

D = 10, E^c^ = 2, E^^<i + E^^d = 9: 

{W^o, W^}, {W2, W,}, {2Wo, Kf}, 

{Wo, W^, Xe}, {2W2, K^}, {Wo, W^,K^}, 

{W2, W^,K^}, {2W^, Kf}, {Wo, W2, Kf, K,}, 

{2Wo, kI\ K,}, {2Wo, K,, K,}, {Wo, W,, xi^^}. 
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{21^0, iff, Ki^\ 4% 
e) D = 12, Y^^d = 0, J2dUd = 11: 

{Kl\Ki^\ K,}. (A.18) 

The indices j, k take admissible values according to Eq. ( |A.13| ). Again, we have carefully 
checked that terms linear in S or S* do not survive when the product of all vertices 
within each set listed in Eqs. ( |A.15| )-( A.18| ) is formed by formally replacing the bilinears 
H\Hi, hIH2 and S*S with 1. 

Nevertheless, at the seven-loop level, we can still construct tadpole supergraphs 
by combining the renormalizable superpotential vertex Wq four times with the higher- 
dimensional Kahler-potential vertices K^^'^ and Kq in Eq. ( |A.13| ). In other words, the set 
of vertices 

{4Wo, Ki'\ K,} (A.19) 

gives rise to the typical seven-loop tadpole graph of Fig. |l](b). Finally, we can check that 
the global constraint of Eq. ( |A.4| ) is satisfied, with n = 7, D = 10 = even, J2d = 4, and 
Ed dVd + Ed dUd = J2ddUd = 9. 



B The Peccei— Quinn-symmetric limit 

In this appendix we shall derive the analytic expressions for the Higgs-boson masses and 
couplings pertinent to the two gauge-singlet SUSY extensions of the MSSM in the PQ- 
symmetric limit, i.e. n/X, ts, 0. Of course, a kinematic situation close to the 

PQ-symmetric limit can more naturally be realized in the NMSSM rather than in the 
MNSSM where Xts/fi is expected to be unsuppressed of order Mgusv Additionally, we 
shall assume that Mh+ ^ M^. For notational simplicity, we have everywhere dropped the 
superscript (0), e.g. from M'f^+\ M^'^^^ etc., as all quantities involved in this appendix are 
evaluated at the tree level. 



In the limit of a heavy H^, the quantity 6 = wA^f^/(2M^) defined in Eq. ([4.14|) is 



much less than 1 and therefore serves as an expansion parameter in our calculations. In 
this limit, it is a reasonable approximation to set fi to its value in the middle of the allowed 
/i^-interval determined by Eq. (|4.13|) , i.e. /x^ = /i^j^ = s^c^Mf, at which Mh-^ is expected 



to approximately acquire its maximum. Then, the tree-level CP-even Higgs-boson mass 
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matrix Mj may be cast, up to terms of order 5 M^, into the approximately diagonal form: 







V 








\ 

i X'v' - Ml) .2,C2, M! + \ XW,, - (i A^.;^ - M|) j 

(B.l) 



by virtue of the orthogonal matrix 



O 



H 



( sign (A/i) 5S/3C2/3 C/3 (1 - I ^^c^^) \ 

-sign (A/i) 5c/3C2/3 S/3 C/3 (1 - i 



V 



1 _ 1 A2p2 
i 2 " '^2/3 



(B.2) 



sign (A/i) 5c2;3 / 

where we have used the short-hand notation S2/3 = sin 2/3 and C2/3 = cos 2/3. 

Likewise, the orthogonal matrix O^, which diagonalizes the CP-odd Higgs-boson mass 
matrix M|, in the PQ-symmetric limit, is easily found to be 

1 / —sign (A/i) 5 1 
VTTP \ -1 -sign (A/i) (5 



:b.3) 



In the PQ-symmetric limit, the CP-odd mass matrix Mp has one massless eigenstate Ai 
and one massive one A2, with = + ^ X^v"^. 

Substituting Eqs. ( [B.2| ) and ( |B.3| ) into Eqs. ( p.51| ) and ( |3.52| ), we obtain, up to order 
(5^, all the couplings HiZZ and HiAjZ: 



9hiZZ 

dHiAiZ 

9HIA2Z 



, 9H2ZZ — 1, QhzZZ 

— S'^ C2/3 , gHiAiZ = , 



9H3AIZ 



sign (A/i) 5 C2/3 , fi-HaAaZ = , gn-^AiZ = 1-^(1 + cjp) 



-sign (A/i) (5 , 
2 



:b.4) 
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Figure 2: Numerical predictions for Mh^ as a function of /i in the MNSSM with = 0, 
for Mh+ = 0.1 (sohd hne), 0.3 (dashed hne), 0.7 (dotted hne), 1 (dash-dotted hne) TeV. 
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Figure 3: Numerical estimates of (a) Mh^ and and of (b) gjj^zz q'h-zZZ^ functions 
of ^ in the MNSSM with ml^ = 0, for Mh+ = 80 (sohd hue), 120 (dashed hue) and 160 
(dotted hue) GeV. 
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Figure 4: The same as in Fig. 0, but with tan j3 = 20. 
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Figure 5: Numerical estimates of (a) Mh^ and and of (b) gjj^zz ^^^^ g'u.^zzi functions 
of /X in the MNSSM with ml^ = 0, for Mh+ = 0.2 (sohd hne), 0.4 (dashed hue), 0.6 (dotted 
hne) and 0.8 (dash-dotted hne) TeV. 
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Figure 6: Numerical predictions for (a) Mh^ and and for (b) g\^zz ^^^^ g'jj^zzi 

functions of /x in the MNSSM with = 0) for Mh+ = 0.3 (soHd Hne), 0.5 (dashed hne), 
0.7 (dotted hne) and 1 (dash-dotted hne) TeV. 
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Figure 7: Numerical predictions for (a) Mh^-, and M^^, and for (b) g\^zz ^^^^ 9h2ZZ^ 
as functions of in the MNSSM, for = 80 (solid line), 120 (dashed line) and 160 

(dotted line) GeV. Numerical estimates of Ma^ are indicated by dash-dotted lines. 
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Figure 8: The maximal predicted value of Mh^ as a function of the charged Higgs-boson 
mass Mh+ in the MNSSM with ml^ = 0. 
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Figure 9: Numerical estimates of (a) Mh^ and and of (b) g'jj^zz ^^"^ 9H2ZZ functions 
of /X in the NMSSM, for Mii+ = 0.3 (solid line), 0.5 (dashed line), 0.7 (dotted line) and 1 
(dash-dotted line) TeV. 
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Figure 10: The same as in Fig. ^ but with k = 0.1. 
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Figure 11: The same as in Fig. |, but with k = -0.1 and KfiA^/X = -0.0021 TeV^. 
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Figure 12: Numerical estimates of (a) Mh^ and and of (b) g'jj^zz ^"^^ Su^zz functions 
of ^ in the NMSSM, with the constraint Ma^ 0, for M^h+ = 120 (sohd hne), 400 (dashed 
hue) and 800 (dotted hne) GeV. 
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